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SUMMARY 


The objective of this study was the evaluation, through both experimenta- 
tion and a literature review, of several advanced concepts for storing natural 
gas at reduced pressure. The advanced concepts included in this study were 
adsorption on high surface area carbon, adsorption in high porosity zeolite, 
storage in clathration compounds, and storage by dissolution in liquid sol- 
vents. The literature review indicated that high storage capacity could be 
obtained with adsorption systems. 

Seventeen carbon samples and seven zeolite samples were then secured and 
evaluated experimentally in a pressurized microbalance apparatus to determine 
their methane adsorption isotherm. The methane storage capacity of each sam- 
ple was also determined using a bench-scale storage system. Results indicated 
that high surface area carbons with high packing density were the best low 
pressure storage mediums. 

A simple mathematical model was used to compare adsorption storage on a 
state-of-the-art carbon with compression storage. The model indicated that a 
vehicle using adsorption storage of natural gas at 3.6 MPa (500 psig) will 
have 36% of the range, on the EPA city cycle, of a vehicle operating on a com- 
pression storage system having the same physical size and a peak storage pres- 
sure of 21 MPa (3000 psig). However, preliminary experiments and current 
literature suggest that the storage capacity of state-of-the-art carbons could 
be improved by as much as 50%, and that adsorption systems having a capacity 
equal to compression storage at 14 MPa are possible. 
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INTRODUCTION 


Recently there has been increased interest in the use of natural gas 
(methane) as fuel for vehicles in the United States, The primary reason for 
this interest is that natural gas is currently lower priced than either gaso- 
line or diesel fuel, but other factors are also important* For example, 
methane-fueled vehicles would decrease both oil imports and vehicular contri- 
bution to photochemical smog. At the present time, there are nearly 30,000 
natural gas fueled vehicles (NGV) in regular operation in the United States; 
most of these are part of commercial fleet operations* However, this number 
pales in significance when compared to the approximately 140,000,000 licensed 
vehicles in the U*S. To achieve a more wide spread acceptance of NGVs and a 
deeper market penetration, more economical methods must be developed for re- 
fueling and for on-board storage of natural gas. 

Two alternative storage approaches are available to provide compressed 
natural gas (CNG)-fueled vehicles with driving ranges comparable to those of 
liquid-fueled vehicles. These are high storage pressure at pressures up to 21 
MPa (3000 psi) and low-pressure adsorption storage at pressures below 3.6 MPa. 
Most NGVs in operation today use high-pressure storage, and many studies are 
available to justify their use — primarily for multi-vehicle fleets. But 
only a limited amount of work has been done in low-pressure absorption storage 
where the biggest potential cost benefits lie. For example, a recent (October 
1982) estimate of costs for converting a 70 vehicle fleet to dual fuel high 
pressure (21 MPa) operation broke down as follows: $83,500 for vehicle con- 

version equipment and $140,000 for a quick fill refueling station. When an 
overnight fill was substituted for the quick fill, refueling station costs 
dropped to $94,000, still more than half the cost of conversion. These costs 
average $3,200 per vehicle for the quick fill option and $2,500 per vehicle 
for an overnight fill refueling station. Payback periods as short as 2 years 
are possible for this fleet if the fleet averages 45,000 miles per vehicle per 
year and the combined average fuel economy is under 20 miles per gallon. 
Clearly, the high cost of compressor equipment is a major stumbling block to 

deeper market penetration. 

• 

The, high costs of refueling stations has also helped to bring about a 
“chicken or the egg" syndrome in this fledgling industry. The private sector 
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will not commit the large capital expenditures required to build public re- 
fueling stations until a substantial market exists® The individual consumer, 
however, will not opt for a dual fuel conversion or a dedicated NGV until 
there is widespread availability of fuel® However, if adequate low pressure, 
on-board storage capacity can be achieved, substantial reductions in the cost 
of a refueling station would result® In fact, a limited number of sites would 
be available, at the outskirts of major metropolitan areas, where vehicles 
could be refueled directly from high pressure (3.6 MPa and above) gas trans- 
mission lines, thus entirely eliminating the need for compressors® Clearly, 
low pressure storage warrants a closer look. 

There have been some efforts in the recent past to store methane at 1*5 

to 3.6 MPa using physical adsorption on solid materials such as carbons and 

1 2 

zeolites. The most recent work, * conducted by Ford Motor Company, has 
emphasized carbon over zeolites. Ford determined that a Union Carbide carbon 
designated 9LXC was the best adsorbent of those investigated. 

The purpose of this present effort was to conduct fundamental studies on 
carbons with high surface areas and zeolites with high internal porosities. 
These studies identified the effect on methane storage of critical sorbent 
characteristics such as surface area, pore size, pore size distribution, 
particle size, and sorbent packing density. 

The work in this program was divided into three major tasks. Task 1 
evaluates the effect of surface area, pore size, pore size distribution, 
particle size, and packing density on the methane storage capacity of carbons 
and zeolites. Task 2 completed a literature review that was already underway 
at IGT and addressed alternate low pressure storage concepts utilizing clath- 
ration compounds and solvents. Task 3 summarized the entire program and 
recommends future research and development. 
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TASK 1. ADSORPTION STUDIES 


Task 1.1 Effect of Storage Medium Characteristics on Methane Storage Capacity 

The objective of this task was to document the effect of surface area, 
pore volume, pore size distribution, particle size, and packing density on the 
methane storage capacity of carbons and zeolites at room temperature and at 
pressures up to 3.6 MPa. 

A. Experimental Efforts — Carbons 
1 . Suppliers 

Samples of high-surface-area carbon blacks and activated carbons were 
obtained from various manufacturers. Table 1 is a list of suppliers and 
product names. 


Table 1. LIST OF CARBONS BY SUPPLIER 


Supplier 

Gulf Oil Chemicals Co. 
(Shawinigan Products Div.) 

Cabot, Corp. 

ICI Americas, Inc. 

Westvaco 

J. T. Baker Co. 

CECA, Inc. 

Witco Chemicals 
Calgon, Corp. 

North American Carbon 


Sample 


Designation 

No, 

Acetylene Black 
50% Dense 

Cl 

100% Dense 

C2 

CSX-179-B 

C3 

DXL-0-8334 

C4 

DARCO SG 

C5 

Nuchar WV-B 

C6 

Nuchar WV-G 

C7 

Acid-Washed Carbon 

C8 

GAC 50G 

C9 

JXC 4 x 6 

CIO 

PCB 30 x 140 

Cll 

BPL 30 x 140 

C12 

G101 

C13 

G104 

C14 

G210 

C15 

G216 

C16 

9LXC 

C17 


Union Carbide 

(No longer in production) 
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2 • Residue After Ignition 

One-gram samples of each of these 17 carbons were placed in weighed, 
porcelain crucibles and heated overnight in air at 1000°C. Weight loss was 
determined, and where possible the residue present after ignition was re- 
covered. Some of the activated carbons obtained commercially are in the form 
of cylindrical pellets and contain a clay binder* Others have been impreg- 
nated with inorganic materials to enhance their ability to absorb reactive 
gases such as hydrogen cyanide. Still others, because of their origin, con- 
tain alkali metal carbonates or alkaline earth phosphates and sulfates. It is 
important to know how much of each sample is carbon in order that only high 
purity carbons might be used to evaluate the importance of surface area and 
pore volume. Such information will be essential in the modeling effort to be 
undertaken in later stages of Task 1. Results are listed in Table 2. 

3 • Acidity 

The key to the adsorption of gases on high-surface-area activated carbons 
is the nature of the carbon surface. The activation of charcoal, for example, 
by partial oxidation in air produces an unstable surface layer. Upon contact 
with air at room temperature, this unstable layer reverts to a more stable 
modification called oxidized charcoal. This oxidized charcoal possesses a 
variety of carboxyl and carbonyl groups on its surface that contribute sig- 
nificantly to its adsorption properties. Boem, for example, has characterized 
the surface groups on oxidized charcoal by titration with bases of different 

o 

strengths and by chemical analysis. He proposes that oxidized charcoal has 
at least four types of surface structures (Figure 1). Structures II and III 
both possess carboxylic acid groups while Structure IV has positively charged 
counter ions associated with it. When mixed with water, Structures II and III 
can dissociate to give a distinctly acid character to the water; Structure IV 
can hydrolyze to yield a base. Although quantitative chemical characteriza- 
tion of the 17 carbon samples with respect to the nature of the surface was 
beyond the scope of the present program, some indication of the surface func- 
tionality was obtained by mixing the carbon with water and measuring the pH. 

One-gram samples of each of the 17 carbons were placed in approximately 
50 mL of degassed, deionized water and mixed in a high-speed blender for 1 
minute under a nitrogen atmosphere. The resulting suspension was allowed to 
settle for 10 minutes, and the pH of the solution was then determined using a 
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Table 2. ANALYTICAL DATA FOR CARBON SAMPLES 


Sample No, 

Residue after 
Ignition, % 

Residue Color 

pH of Carbon 
Sample 

Cl 

<0.5 

— 

5.3 

C2 

<0.5 

— 

5.4 

C3 

<0.3 

— 

10.1 

C4 

7.6 

Orange 

8.4 

C5 

14.5 

Cream 

8.1 

C6 

3.9 . 

Colorless Glass 

3.8 

Cl 

6.6 

Orange 

10.0 

C8 

<0.3 

— 

6.3 

C9 

13.5 

White 

8.3 

CIO 

<0.1 

— 

9.8 

Cll 

82.0 

Orange 

9.6 

C12 

73.6 

Colorless Glass 

9.7 

C13 

11.8 

Colorless Glass 

6.4 

C14 

12.0 

Colorless Glass 

2.6 

C15 

1.5 

Colorless Glass 

9.9 

C16 

1.2 

Colorless Glass 

9.6 

. C17 

1.0 

Black Fibers 

9.9 
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m 12 

Figure 1. FOUR TYPES OF SURFACE STRUCTURES FOUND IN OXIDIZED 
CHARCOAL AS PROPOSED BY BOEM IN REFERENCE 3 
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Corning Model 476223 Semi Micro Combined Electrode and a Markson Scientific 
Model 90 digital pH meter* The results are listed in Table 2* It should be 
noted that with those carbons having a high amount of residue after ignition, 
a substantial contribution to pH could come from sources other than carbon* 


4* Surface Area Determination — Theory 


The physical adsorption of gas molecules onto the surface of a solid 
results from relatively weak interaction between molecules of the solid and 
molecules of the gas* This interaction is frequently referred to as van der 
Waal’s forces* Physical adsorption contrasts with chemical adsorption, where 
chemical bonds are responsible for the interaction between solid and gas* In 
physical adsorption, the quantity of adsorbed gas increases with decreasing 
temperature. In chemical adsorption, because it depends upon chemical reac- 
tion, the quantity of adsorbed gas decreases with decreasing temperature. For 
this reason, adsorption measurements to determine surface area or pore distri- 
bution are made at lower temperatures where physical adsorption predominates* 
Brunauer, Emmett, and Teller (B.E.T.)^* 5 have derived the following expression 
for relating the volume of gas adsorbed to the nature of the solid surface and 
the temperature and pressure of the gas: 


V CP 

y = 2 

a (P o - P) [1 + (C - 1) P/P J 
s s 

where — 


V. 

P 

P, 


m 


Volume of adsorbed gas 

Volume of gas required to form a monolayer over entire surface 
Pressure of gas 

Equilibrium vapor pressure of gas and its liquid at temperature of 
measurement (760 torr at 77.35°K for nitrogen) 

A constant dependent upon the nature of the solid adsorbent. 


This expression describes the great majority of low-temperature adsorption 
data. Vm can be calculated from a series of physical measurements of the 
volume of gas adsorbed as a function of pressure at a fixed temperature by 
rearranging Equation 1 to a linear form: 
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( 2 ) 


P 1 C - 1 P 

V(P-P) VCVC P 
as m m s 

A plot of P/V fl (P g - P) versus P/P g gives a straight line with an intercept at 
l/V m C and a slope of (C - l)/V m C. The value of V m is readily determined from 
a series of such measurements. The linearity of the B.E.T. equation can be 
expected to hold only in the region of pressures where P/P_ is between 0.05 
and 0.3. A more detailed discussion of B.E.T. theory can be found in Refer- 
ences 4 and 5. 

5. Surface Area Determination — Practice 

The actual surface area was determined on an ORR Surface-Area Pore-Volume 
Analyzer. The procedure used is as follows. 

A weighed carbon sample is placed in a small containment vessel and 
attached through a valve to a gas manifold of known volume, V^. The sample is 
warmed to over 100°C under a vacuum (<10 microns), then sealed off from the 
manifold and cooled to room temperature. The manifold is charged with helium 
to a known pressure, P^. The containment vessel with the sample is then 
opened to the manifold, where upon the manifold pressure drops to P£« The 
volume of the containment vessel with sample, V^, can now be calculated as 
follows: 


V 


C 


P V 
* 1 V 1 


- V, 


(3) 


This procedure assumes that no helium adsorbs on the sample at these tempera- 
tures. The containment vessel and sample are again evacuated and then cooled 
with liquid nitrogen. In a manner similar to that above, the manifold is 
pressurized again but with nitrogen at room temperature, T R , and the pressure, 
P^, is accurately determined. The manifold is then opened to the containment 
vessel, allowed to come to equilibrium, and the new pressure, P 2 , is measured. 
Using Equation 4, the volume of gas adsorbed, V a , at pressure P£ can be calcu- 
lated. 


P V 
*1 V 1 


P V P V 
Vl , 2 v c 760 

T 773k 273 a 


(4) 
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The process is then repeated for several incremental additions of nitrogen 
gas* with additional terms being added to Equation 4 to account for the fact 
that the initial pressure is no longer zero. In this way* several simultane- 
ous values for V a and P 2 can be generated from which a B.E.T. plot can be con- 
structed using Equation 2. 

The above discussion is somewhat simplified in that we have neglected 
several corrections to V a that must be included because of the nonideal 
behavior of nitrogen and because of peculiarities of the instrument itself. 

But once a corrected B.E.T. plot is made using the various V a versus ^^s 
data points, the specific surface area of the sample may be calculated from 
the following formula: 

Specific Surface Area - tj.v. (Slope' ^Intercept) (5) 


where — 

7 —90 7 

M.A. = Molecular area in era (16.2 X 10~^ u cm for nitrogen) 
N = Avogadro’s number 

M.V. = Molar volume in cm 

This equation can be simplified to Equation 6 for nitrogen: 


Specific Surface Area = 


4.35 

Slope + Intercept 


( 6 ) 


The results for the specific surface area determinations for the 17 carbon 
samples are listed in Table 3. It should be noted that with high surface area 
samples where V m is large, the slope and the intercept become small, thus in- 
creasing measurement error. Although the specific surface areas are given to 
three significant figures in Table 3, the relative error for these measure- 
ments may be as much as ±10%. 

6. Packing Density 

The packing density of the adsorbent, defined here as the mass in grams 
of one cubic centimeter of settled material, is one of the critical parameters 
associated with the use of adsorbents for CNG storage in automotive vehicles. 
Packing density can be conveniently measured by placing a weighed amount of 
material in a graduated cylinder and vibrating the cylinder a rate of 100 Hz 

10 



Table 3. SPECIFIC SURFACE AREA, APPARENT DENSITY, AND 
SURFACE AREA PER LITER FOR 17 CARBON SAMPLES 


Sample 

Designation 

Specific 

Surface Area, nr/g 

Packing 
Density, g/ cnr 

Surface Area 
km^/L 

Cl 

76 

0.10 

0.008 

C2 

74 

0.20 

0.015 

C3 

1600 

0.13 

0.21 

C4 

1030 

0.44 

0.45 

C5 

700 

0.45 

0.32 

C6 

1610 

0.30 

0.48 

Cl 

1260 

0.45 

0.57 

C8 

480 

0.37 

0.18 

C9 

1030 

0.56 

0.58 

CIO 

1050 

0.45 

0.47 

Cll 

1270 

0.44 

0.56 

C12 

1100 

0.47 

0.52 

C13 

1680 

0.30 

0.50 

C14 

1650 

0.30 

0.50 

C15 

1420 

0.50 

0.71 

C16 

1370 

0.50 . 

0.69 

C17 

1280 

0.32 

0.41 


11 



or higher for several minutes to allow the sample to settle under the influ- 
ence of gravity. Packing density measured in this way can be used to calcu- 
late the quantity of material that can be placed in a storage cylinder from a 
free-flowing reservoir. Table 3 lists the packing density of the 17 carbon 
samples in grams per cubic centimeter. The specific surface area and the 
packing density can also be used to estimate the surface area available for 
storage in a 1-liter container filled with carbon. The final column of 
Table 3 contains such estimates for the various carbon samples. 

7 . Particle Size — Activated Carbon Pellets 

All of the activated carbons supplied to us were in the form of pellets 
or granules while the carbon blacks were in the form of powders. In the case 
of the activated carbons, the particle sizes fell within narrow limits which 
is likely the result of a screening step in their manufacture. The various 
suppliers provided us with information on the upper and lower limits of par- 
ticle size for each sample as measured in U.S. mesh. Table 4 is a particle 
size conversion table comparing the U.S. mesh size with the equivalent par- 
ticle diameter in microns. This table is included solely for convenience. 
Table 5 contains a list of the 17 carbon samples and the maximum and minimum 
mesh sizes for the 13 activated carbons. Because better than 95% of the par- 
ticles in any given sample fell within the narrow limits of size listed by the 
respective supplier, no further attempt was made to determine a particle size 
distribution for the activated carbon samples. 

8. Particle Size Distribution — Carbon Black Powders 

The particle size distributions for the four carbon black samples were 
determined on an automated Coulter counter. Table 6 lists the results for 
50% dense Shawinigan Acetylene Black. These results do not agree with the 
particle size distribution provided by the manufacturer. When viewed under an 
electron microscope, Shawinigan Acetylene Black appears to consist of spheri- 
cal particles clumped together in much larger agglomerates. The median par- 
ticle size for these spheres as determined from a micrograph is 42.5 nm 
whereas the mean agglomerate size as determined by the Coulter counter is 23 
urn or about 540 times larger. Since it is the agglomerate size rather than 
the particle size which determines handling characteristics such as packing 
density, the agglomerate size distribution is the more important parameter. 
Table 7 lists the agglomerate size distribution for 100% dense Shawinigan 
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Table 4. PARTICLE SIZE CONVERSION TABLE 


Mesh Approximate Size 

Size in Microns* 


4 

4760 

6 

3360 

8 

2380 

12 

1680 

16 

1190 

20 

840 

30 

590 

40 

420 

50 

297 

60 

250 

70 

210 

80 

177 

100 

149 

140 

105 

200 

74 

230 

62 

270 

53 

325 

44 

400 

37 

625 

20 

1250 

10 

2500 

5 


* 1 mm = 1000 microns. 
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Table 5. PARTICLE SIZE DISTRIBUTION BY MESH SIZE 


Number 

Carbon 


Size 

Cl 

Shawinigan 50% 

(see 

Table 6) 

C2 

Shawinigan 100% 

(see 

Table 7) 

C3 

Cabot CSX-179-B 

(see 

Table 8) 

C4 

ICI DXL-0-8334 

20 x 

60 mesh 

C5 

ICI DARCO-SG 

20 x 

60 mesh 

C6 

Nuchar WV-B 

14 x 

35 mesh 

C7 

Nuchar WV-G 

12 x 

14 mesh 

C8 

Baker Acid Washed 

(see 

Table 9) 

C9 

CECA GAC-50G 

20 x 

50 mesh 

CIO 

Witcarb JXC 

4 x 

6 mesh 

Cl 1 

Calgon PCB 

30 x 

140 mesh 

C12 

Calgon BPL 

30 x 

140 mesh 

C13 

North American G101 

10 x 

25 mesh 

C14 

North American G104 

14 x 

35 mesh 

C15 

North American G210 

8 x 

16 mesh 

C16 

North American G216 

14 x 

35 mesh 

C17 

Union Carbide 9LXC 

12 x 

28 mesh 
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Table 6. PARTICLE SIZE DISTRIBUTION FOR 50% DENSE 
SHAWINIGAN ACETYLENE BLACK 



Less than 4 microns 

4.0 to 5.0 microns 

5.0 to 6.4 microns 
6.4 to 8.0 microns 

8.0 to 10.1 microns 

10.1 to 12.7 microns 
12.7 to 16.0 microns 

16.0 to 20.2 microns 

20.2 to 25.4 microns 
25.4 to 32.0 microns 

32.0 to 40.3 microns 

40.3 to 50.8 microns 
Larger than 50.8 microns 



0 

1.1 

2.4 

4.4 

5.9 

6.5 

7.6 
11.2 
16.0 
20.4 
13.7 

6.9 

3.9 


15 



Table 7. PARTICLE SIZE DISTRIBUTION OF 100% DENSE 
SHAWINIGAN ACETYLENE BLACK 


Size Range 

Less than 4 microtis 

4.0 to 5.0 microns 

5.0 to 6.4 microns 
6.4 to 8.0 microns 

8.0 to 10.1 microns 

10.1 to 12.7 microns 
12.7 to 16.0 microns 

16.0 to 20.2 microns 

20.2 to 25.4 microns 
25.4 to 32.0 microns 

32.0 to 40.3 microns 

40.3 to 50.8 microns 
Larger than 50.8 microns 


Weight Percent 
0 

5.9 
10.3 

10.3 

9.9 

9.8 
10.1 
11.2 

14.4 
10.7 

4.8 
1.6 
0 
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Acetylene Black; Table 8, for Cabot Corporation carbon black CSX-179-B; and 
Table 9, the agglomerate size distribution for Baker-Acid Washed Carbon* 

9# Cumulative Pore Volume 

Both nitrogen adsorption and desorption isotherms were determined at 78°K 
for all carbon samples and the results used to estimate the cumulative pore 
volume of each sample and to ascertain the pore size distribution. Table 10 
lists 16 carbon samples and the cumulative pore volume of all pores less than 
200 Angstroms in radius. It should be noted that our cumulative numbers have 
a cut-off for maximum pore radius at 200 Angstroms and, therefore, pores larger 
than 200 A are ignored, since the interior surfaces of pores larger than 200 A 
in radius would be little different from flat surfaces. However, as a conse- 
quence, our cumulative pore volume may differ somewhat from the total pore 
volume numbers advertised for these carbons by their various manufacturers. 

It can be seen from Table 10 that Cabot f s CSX-179-B has the highest cumu- 
lative pore volume of the carbon blacks, while Nuchar WV-B is highest among 
the activated carbons. Of particular interest are the four North American 
Carbon, Inc., samples. The 100 series samples have high cumulative pore 
volume while the 200 series have the lowest measured. 

10. Pore Size Distribution 

The nitrogen desorption isotherm data were used to plot pore size distri- 
butions for 16 of the carbon samples. Subsequent work indicated no identifi- 
able relationship between pore size distribution and methane storage capacity. 
Individual plots of pore size distribution as pore volume versus pore radius 
are provided in Appendix A for carbons Cl through C16. 

11. Methane Adsorption Isotherms for Carbons 

Figure 2 is a schematic diagram of the apparatus used to determine the 
methane adsorption isotherms for selected carbon and zeolite samples. The 
experimental technique used was to suspend a measured quantity of adsorbent 
from an electronic transducer type balance in a methane atmosphere and deter- 
mine the weight of the sample as the vessel is pressurized from 101.4 KPa (0 
psig) to 3.6 MPa (500 psig) in 50 psi steps. Weight changes as little as 1 
part in 2000 can be measured with this apparatus. However, it is necessary to 
make corrections on the indicated weight for the buoyancy of the sample. The 
following may serve as an example. In the carbon experiments, a sample of JXC 
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Table 8. PARTICLE SIZE DISTRIBUTION FOR 
CABOT CSX-179-B 


Size Range 


Weight Percent 


Less 

than 4 microns 

0 

4.0 

to 

5.0 

microns 

7.7 

5.0 

to 

6.4 

microns 

12.5 

6.4 

to 

8.0 microns 

14.0 

8.0 

to 

10.1 

microns 

15.4 

10.1 

to 

12.7 

microns 

14.7 

12.7 

to 

16.0 

microns 

11.8 

16.0 

to 

20.2 

microns 

8.2 

20.2 

to 

25.4 

microns 

4.9 

25.4 

to 

32.0 

microns 

2.7 

32.0 

to 

40.3 

microns 

1.5 

40.3 

to 

50.8 

microns 

0.8 


Larger than 50.8 microns 5.8 
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Table 9. PARTICLE SIZE DISTRIBUTION FOR BAKER ACID WASHED CARBON 
Size Range Weight Percent 


Less 

than 3.2 

microns 

0 

3.2 

to 

4.0 

microns 

1.3 

4.0 

to 

5.0 

microns 

2.7 

5.0 

to 

6.4 

microns 

3.4 

6.4 

to 

8.0 

microns 

4.4 

8.0 

to 

10.1 

microns 

4.2 

10.1 

to 

12.7 

microns 

2.4 

12.7 

to 

16.0 

microns 

1.3 

16.0 

to 

20.2 

microns 

1.1 

20.2 

to 

25.4 

microns 

0.7 

25.4 

to 

32.0 

microns 

0.5 

32.0 

to 

40.3 

microns 

0.2 

40.3 

to 

75 

microns 

17.5 

75 

to 

100 

microns 

35.3 

150 

to 

300 

microns 

17.7 

300 

to 

600 

microns 

6.9 


Larger than 600 microns 
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Table 10. 'CUMULATIVE PORE VOLUME FOR PORES OF 
RADIUS LESS TRAN 200 ANGSTROMS 


Number 

Carbon 

Pore Volume 

Cl 

Shawinigan 50% 

0.10 ml/g 

C2 

Shawinigan 100% 

0.11 ml/g 

C3 

Cabot CSC-179-B 

0.99 ml/g 

C4 

ICI DXL-0-8334 

0.17 ml/g 

C5 

ICI DARCO-SG 

0.53 ml/g 

C6 

Nuchar WV-B 

0.66 ml/g 

Cl 

Nuchar WV-G 

0.14 ml/g 

C8 

Baker Acid Washed 

0.82 ml/g 

C9 

CECA GAC-50G 

0.09 ml/g 

CIO 

Witcarb JXC 

0.13 ml/g 

Cll 

Calgon P.CB 

0.05 ml/g 

C12 

Calgon BPL 

0.09 ml/g 

C13 

North American G101 

0.54 ml/g 

C14 

North American G104 

0.58 ml/g 

C15 

North American G210 

0.024 ml/g 

C16 

North American G216 

0.026 ml/g 
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Figure 2. EXPERIMENTAL APPARATUS USED TO DETERMINE 
METHANE ADSORPTION ISOTHERM 
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carbon weighing 10*076 grains under one atmosphere pressure (101*3 KPa) of 
methane was pressurized to 35 atmospheres (3.6 MPa). Using the adsorption 
isotherm data previously obtained, we would expect this sample to gain 0*856 
gram in weight due to methane adsorption* The microbalance, however, indi- 
cated a gain of only 0.746 gram. The difference, 0.11 gram, is due to 
buoyancy effects. The carbon in a 10.076 gram sample of JXC carbon has an 
actual volume of only 4.80 cc and displaces 0.110 gram of methane at 35 atmo- 
spheres. The key to calculating the buoyancy correction is the actual density 
of the solid material. Witco’s JXC carbon (CIO) has a packing density of only 
450 grams/liter; however, almost 80% of a sample of JXC is accessible void 
space which can fill with compressed gas and does not contribute to buoyancy. 
In the case of most of the carbon samples, the solid material contributing to 
buoyancy effects can be considered to have an effective density equal to that 
of graphite, 2100 grams/liter. 

Adsorption isotherms for each of the carbons C3 through Cl 7 are provided 
in Appendix B. Table 11 summarizes the results. Since Ford Motor Co. 1 and 
others are working with adsorption systems operating at about 2.2 MPa (300 
psig), we have reported data for methane adsorption at this pressure in Table 
11. We have also included data for methane adsorption at 3.6 MPa (500 psig) 
which is the maximum pressure used in the present study. Because of the par- 
ticulate nature of the Shawinagan Acetylene Blacks, these carbons tended to 
fluidize during depressurization of the microbalance and reliable adsorption 
measurements were not possible for Samples Cl and C2. 

B. Experimental Effort — Molecular Sieves 
1 • Description of Samples 

Zeolites . Six samples of commercially available crystalline molecular 
sieves were obtained from the Linde Division of Union Carbide,. Four of these 
samples are zeolites of the Sodalite Group, one sample is of the Chabazite 
Group, and one sample is of the Mordenite Group. One additional sample was 
prepared by ion exchange. Table 12 lists the samples by experimental 
designation (Zl through Z7) and provides information on structure, critical 
diameter, and nonframework cation. 
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Table 11. ADSORPTION CHARACTERISTICS OF CARBON SAMPLES* 


Sample 

Packing 

No. 

Density 

C3 

0.15 

C4 

0.42 

C5 

0.44 

C6 

0.26 

Cl 

0.49 

C8 

0.40 

C9 

0.54 

CIO 

0.45 

Cll 

0.44 

C12 

0.47 

C13 

0.24 

C14 

0.27 

C15 

0.49 

C16 

0.45 

C17 

0.32 


* At 23°C. 


Methane Adsorbed in 


;rams/gram 


2.2 MPa 3.6 MPa 


0.099 

0.110 

0.074 

0.081 

0.050 

0.053 

0.072 

0.085 

0.085 

0.091 

0.042 

0.045 

0.069 

0.076 

0.083 

0.087 

0.090 

0.095 

0.074 

0.079 

0.079 

0.091 

0.068 

0.078 

0.092 

0.099 

0.095 

0.103 

0.101 

0.115 
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Table 12. INFORMATION ON VARIOUS ZEOLITES 

Linde Critical Nonframework 


Designation 

Structure 

Designation 

Diameter, A 

Cation 

Z1 

Zeolite A 

3A 

3 

Potassium 

Z2 

Zeolite A 

4A 

4 

Sodium 

Z3 

Zeolite A 

5A 

5 

Calcium 

Z4 

Zeolite X 

13X 

10 

Sodium 

Z5 

Zeolite A 

— 

n 

Lithium 

Z6 

^lordenite 

AW-300 

3 to 4 

Mixed 

Z7 

Chabazite 

AW- 5 00 

4 to 5 

Mixed 


The structure of the Sodalite Group, Zeolite A, is based on frameworks 
that are simple arrangements of truncated octahedra. These truncated octa- 
hedra share square and hexagonal faces. In the structure of Zeolite A, the 
octahedra are linked by adjoining cubes (Figure 3). This produces a central 
truncated octahedron with an internal cavity of 11 A diameter. Access to this 
cavity is by way of the six apertures, which are the cubes with a free 
diameter of 4.2 A. In Zeolite A of the formula Na£ [ (AK> 2 2 (^ 02 ) 12 ^ *27 H 2 O, 
there are 12 nonframework sodium ions per unit cell. Eight of these sodium 
ions reside in the center of the eight hexagonal faces and are referred to as 
Type I cations. The other four cations, Type II, occupy positions adjacent to 
the openings that interconnect the cavities. When completely hydrated, these 
four ions probably float within a coordination sphere of water molecules; but 
when dehydrated, they locate on the walls of the cavity and exert an influence 
on the critical diameter of the opening between cavities. The critical 
diameter is, for our purposes, the largest diameter a gas molecule can have 
and still pass between cavities. As is illustrated in Table 12, the large 
diameter of the potassium ion (2. 66 A) restricts the critical diameter in Type 
A zeolites to 3A. Thus, methane with a kinetic diameter of 3.8A cannot gain 
access to the interior of the cavity, whereas water molecules with a kinetic 
diameter of 2.65A can be admitted. Type A zeolites with smaller nonframework 
ions (sodium and lithium, having ionic diameters of 1.90A and 1.36A, respec- 
tively) increase the critical diameter to 4A and thus permit methane to enter 
the cavity. Dipositive calcium ions increase the critical diameter to 5 A 
partly because of size, 1.98A, and partly by perturbing the framework itself. 
The Type X zeolites, on the other hand, have a different crystallographic 
morphology in which adjacent cavities share large hexagonal openings. Type X 
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Figure 3. CUBIC ARRAY OF TRUNCATED OCTAHEDRA 
IN TYPE A ZEOLITE 
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zeolites have a critical diameter of 10A. In dehydrated chabazite, the cavi- 
ties share ellipsoidal apertures with dimensions of 4.4A by 3.1A. The zeolite 
mordenite possesses a predominance of 5-member rings that align to form circu- 
lar channels with a free diameter of 6.6A. This structure is sensitive to 
stacking faults, however, which can reduce the large-channel diameter from 6,6 
to 4A. A more detailed discussion of crystalline molecular sieves can be 
found in Reference 6, Of importance here is that the nonframework cations in 
both Type A and Type X zeolites can easily be exchanged* Thus, lithium- 
containing Type A zeolite can be prepared from sodium-containing material 
(Linde 4A). 

The phenomenon of ion exchange was used to prepare Sample Z5* The fol- 
lowing procedure was used* 

One hundred grams of Linde 4A molecular sieves were placed in a 1-liter 
flask fitted with a reflux condenser, and 500 niL of a 0.1M lithium chloride 
solution was added. The mixture was heated to 80 °C and maintained at this 
temperature for 24 hours with occasional gentle agitation. After 24 hours, 
the LiCl solution was drained and replaced with 500 mL of fresh LiCl solution 
and reheated. This process was repeated through four consecutive 24-hour 
treatments. The zeolite was then recovered and washed 10 times with 500 mL 
each of deionized water, drained, and allowed to air dry for 24 hours at room 
temperature. The lithium-containing zeolite was then dried at 200°C in air 
for 1 week, at 360°C in air for 24 hours, and at 360°C for 24 hours under 
vacuum (<10 microns). The vacuum flask was finally cooled and pressurized 
with dry argon to 1 atm. The lithium-containing zeolite was then loaded into 
an airtight container under an argon atmosphere and stored for later use. 

2. Adsorption Isotherms 

All zeolite samples were dried prior to determining their adsorption 
isotherm. The drying process used was that of heating the sample to 360°C for 
24 hours under a continually increasing vacuum which was below 10 microns 
pressure at the end of the 24 hour period. The samples were cooled under 
vacuum and then pressurized to 1 atmosphere with dry argon. The sample was 
transported to the microbalance under argon and quickly transferred to the 
sample pan in a manner designed to limit its exposure to moisture laden air 
(ambient laboratory conditions) to a few minutes. 
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A problem was encountered when attempts were made to calculate a correc- 
tion factor for buoyancy* 

In the case of the zeolites, calculations are not easy. For example, 
when a sample of activated 13X zeolite pellets is packed into a graduated 
cylinder and weighed, it exhibits a packing density of 0.66 g/cc. When a 
single pellet is weighed and its geometrical volume is determined, an apparent 
density of 1.29 g/cc is calculated* However, Linde claims that 13X pellets 
contain an internal void space of 0.36 cc/g. Thus, while a single pellet of 
13X has an apparent density of 1.29 g/cc and one gram would occupy 0.77 cc, 
0.36 cc of this volume is void space and only 0.41 cc is solid material. 

Thus, for purposes of calculating buoyancy effects, the effective density of 
13X zeolite is 2.44 cc/g. However, the aperture size for 3A zeolite is not 
accessible to methane and cannot fill with compressed gas. Consequently, this 
void volume will contribute to buoyancy effects. In 3A zeolite, the aperture 
size is larger than the kinetic diameter of methane and the 0.30 cc/g void 
volume in 5A zeolite can fill with compressed gas and so does not contribute 
to buoyancy. The effective density of 3A zeolites would, therefore, be 1.35 
g/cc while the effective density, for purposes of calculating buoyancy, for 5A 
zeolites would be 2.2 g/cc. 

To summarize, while the apparent densities of pellets of 3A, 5A, and 13X 
zeolites are 1.33 g/cc, 1.33 g/cc, and 1.29 g/cc, respectively, the effective 
density, on theoretical grounds, for use in calculating buoyancy should be 
1.33 g/cc, 2.2 g/cc, and 2.44 g/cc, respectively. Unfortunately, attempts to 
use these theoretical densities to interpret actual experimental data produce 
less than satisfactory results. 

An empirical solution to this problem was arrived at by careful examina- 
tion of the adsorption data for 13X and 5A zeolite samples. These zeolites 
have large critical diameters which are completely accessible to methane. 
Therefore, a straightforward calculation of buoyancy is possible. The adsorp- 
tion isotherm for 13X zeolite exhibits very little change over the incremental 
pressure change between 3.2 MPa (450 psig) and 3.6 MPa (500 psig). The raw 
raicrobalance data, however, does show an appreciable change in weight due to 
changes in buoyancy. If we assume that the other zeolite samples behave in 
the same manner, that is, that the weight change observed between 450 psig and 
500 psig is due almost entirely to changes in buoyancy, we can then use the 
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raw microbalance data as a direct measurement of the buoyancy and eliminate 
the need for a calculated value* 

The relative errors introduced by this assumption are cumulative. In 
other words, a 1% relative error in estimating the buoyancy between 450 psig 
and 500 psig becomes a 10% relative error over the entire range of pressures 
measured (0 to 500 psig)* As a consequence, we cannot have a high degree of 
confidence in the adsorption isotherms for 3A, 4A, AW-300, and AW-400. 

However, since the performance of these materials falls far short of the 
performance of the 5A zeolites, they are effectively eliminated from 
consideration as an advanced concept for methane storage anyway and no harm is 
done* 


Individual adsorption isotherms for the seven zeolite samples are pro- 
vided in Appendix* B* Table 13 is a summary of their performance. 


Table 13. ADSORPTION CHARACTERISTICS OF ZEOLITE SAMPLES 


Sample 

No. 

Packing 
Density 
g/cc _ 

Methane Adsorbed in grams/gram 
2.2 MPa 3.6 MPa 

Z1 

0.77 

0.024 

0.025 

Z2 

0.76 

0.037 

0.039 

Z3 

0.74 

0.047 

0.049 

Z4 

0.66 

0.048 

0.050 

Z5 

0.72 

0.048 

0.051 

Z6 

0.95 

0.038 

0.039 

Z7 

0.79 

0.030 

0.031 

C. Bench- 

■Scale Storage Experiments 




Methane adsorption is only one of several important parameters contrib- 
uting to the performance of an adsorption storage system. Packing density and 
internal pore volume are others. To assess the performance of adsorbents 
under realistic conditions, samples of adsorbent were packed into a stainless 
steel cylinder having an internal volume of 75 cm and the cylinder evacuated 
with a vacuum pump. The evacuated cylinder and sample were weighed on a top 
loading balance accurate to ±0.01 gram. The cylinder was then pressurized to 
3.6 MPa (500 psig) with dry methane. The heat of adsorption was dissipated by 
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blowing air at 25 °C over the cylinder for two hours while maintaining a con- 
stant pressure with a high pressure regulator. The cylinder was then reweighed. 
This provided a direct measure of the total mass of methane stored at 3.6 MPa 
and includes both that portion stored by adsorption and that portion stored by 
compression in the void spaces. The mass (g) of methane stored per unit 
volume (L) is here defined as methane storage capacity. The dependence of 
methane storage capacity upon packing density and specific adsorption is 
illustrated for selected samples in Figure 4. The solid lines represent a 
locus of points determined by the the theoretical relationship between packing 
density and methane storage for carbons having specific adsorptions of 0.165, 
0.115, 0.100, and 0.085 g/g at 3.6 MPa (35 atm or 500 psig) at 25°C. Actual 
experimentally determined values for carbon C6, CIO, C15, C15C, and C17 are 
also plotted in Figure 4. In Reference 7, Amos Golovoy reports success in 
increasing the packing density of 9LXC by crushing the pellets. The storage 
capacity at 3.6 MPa and packing density reported by Golovoy are plotted in 
Figure 7 as Point Rl. The same procedure was applied by IGT to increase the 
apparent density of Sample Cl 5. A 50 gram sample of this carbon was placed in 
a hydraulic press and crushed at a pressure of 1000 psi (7.3 MPa). The broken 
pellets were then placed in a mortar and fragmented further with a pestle 
using a rocking motion rather than a grinding motion. This produced a large 
number of fine particles with a minimum amount of powder. This carbon was 
then dried at 140°C in a vacuum oven and the methane adsorption isotherm was 
determined. Although the crushed carbon sample now exhibited a packing 
density of 0.58 g/cc, the adsorption isotherm for this carbon was identical to 
that of the original indicating that the adsorption characteristics were not 
changed by the crushing process. This carbon sample was designated as C15C 
and is also plotted in Figure 4. 

Barton j5t^ al.^ have determined the methane adsorption isotherm at 25 °C 
for a carbon which was not included in this study, Amoco GX-32. This carbon 

O 

has the highest surface area of any known carbon, 2500 nr/g (Nitrogen BET) and 

the highest specific adsorption for methane of any carbon reported in the 

literature, 0.165 g/g at 3.6 MPa and 25°C. Point R2 in Figure 4 represents 
the storage capacity of Amoco GX-32 and was placed on the basis of information 

provided by Barton. It is possible that the packing density of this material 

also might be increased by crushing. 
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SPECIFIC ADSORPTION, 
g-CH 4 /g-C 



Figure 4. RELATIONSHIP OF SPECIFIC ADSORPTION AND PACKING DENSITY 
OF CARBON TO METHANE STORAGE CAPACITY AT 3.6 MPa PRESSURE 
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The amount of methane stored by compression alone in an empty cylinder at 
14 MPa (2000 psig) would be about 108 grams/liter. 

As can be seen from Figure 4, the storage capacity of Sample C15C is 
superior to any thus far reported in the literature and is 70% that of com- 
pression storage at 14 MPa. If, by some chemical means, the specific adsorp- 
tion of C15C could be enhanced from 0*100 to 0*120 g/cc, at 3*6 MPa the 
storage capacity would also increase by 20%. Such a goal is not unreasonable 
but no method of attaining it has been identified as yet. 

Caution should be used, however, when comparing an adsorption system with 
a compression storage system on the basis of methane storage capacity alone. 
The amount of natural gas stored by both systems is always greater than the 
amount delivered but because of the shape of the adsorption isotherm, this 
latter quantity is not a linear function of pressure for the adsorption system 
as it is for simple compression storage* The most practical method of deter- 
mining the amount of methane that can be delivered by an adsorption system is 
to construct a bench-scale apparatus. This was accomplished by connecting the 
75 cm stainless steel cylinder filled with the appropriate adsorbent to a gas 
manifold and a wet test meter. To determine a practical delivery capacity for 
the various adsorbents, the cylinder described above was pressurized with 
methane and the methane was then slowly bled through a wet test meter over a 
2-hour period* The quantity of methane delivered by the system in cycling 
from 3.6 MPa (500 psig) to 101*4 KPa (0 psig) was determined directly. An air 
stream at 25°C was directed at the cylinder to maintain temperature during 
desorption. Table 14 lists the quantity of methane delivered in grams per 
liter and the percentage of the total adsorbed which this represents for 16 of 
the carbon samples and two of the best zeolite samples. 

D. Discussion of Results 

As can be seen from Table 14 in the previous section. Carbon Sample C15C 
exhibits the best performance, most methane delivered, of all samples tested. 
However, this sample does not have the highest specific adsorption for methane 
nor the highest micropore volume. It does have one of the highest specific 
surface areas and the highest packing density of all carbons tested. These 
two parameters then become the most desirable characteristics of a good sor- 
bent material for an adsorption storage system. 
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Unfortunately, the methods by which the specific surface area of an 
activated carbon can be increased, by the creation of micropores through 
catalytic oxidation or chemical leaching, also decrease the packing density. 

Table 14. METHANE DELIVERED BY ADSORBENT SYSTEMS 
CYCLED FROM 3.6 MPa TO 101.4 KPa 

Sample Grams of Methane per Fraction of Total 

Designation Liter of Storage in Storage 


C3 

34.8 

91.6%’ 

C4 

45.2 

86.8% 

C5 

36.5 

86.1% 

C6 

41.4 

92.8% 

C7 

54.8 

86.8% 

C8 

31.8 

85.1% 

C9 

51.4 

86.7% 

CIO 

49.2 

84.1% 

Cll 

51.8 

84.9% 

C12 

48.4 

86.9% 

C13 

40.1 

93.3% 

C14 

39.4 

94.3% 

C15 

56.3 

84.0% 

C15C 

62.0 

83.8% 

C16 

55.9 

85.5% 

C17 

51.0 

89.4% 

Z3 

44.1 

84.7% 

Z4 

44.1 

87.6% 


Thus, for example, Amoco GX-32 with a high specific surface area also has a 
low packing density. Therefore, improvements in storage capacity of carbons 
like C15C by increasing surface area will likely reduce packing density. 

On the other hand, none of the carbons used in this study have been 
optimized for methane adsorption. The chemical nature of these activated 
carbons, surface groups, or chemical additives, were optimized for the 
adsorption of active gases such as HCN or ^S. A crude index of a carbon's 
ability to adsorb a chemical species can be determined by dividing its 
specific adsorption in g/g at a fixed pressure by the specific surface area in 
cm /g and multiplying by Avogadro's number over the molecular weight of the 
adsorbed species. The resulting quantity is the coverage in terms of 
molecules adsorbed per cm of surface area. For example. Equation 7 is a 
calculation for the coverage of C17 at 3.6 MPa. 
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0,115 g/g 6.023 x 10 23 molecules /mole „ „„ 1/v 14 . , , 2 ,-, x 

A 9 — x r z " "„ ‘ T me TwTI = 3.38 x 10 molecules/cm (7) 

1280 x 10 4 cm 2 /g 16 g rams /® ole 

Table 15 lists the methane coverage for representative carbon samples. 

It can be seen that the coverage for C15 is about average. If the methane 
coverage of C15 can be increased to equal that of C17 or C8, without altering 
the surface area or the packing density, a 35% improvement in methane adsorp- 
tion and a 25% improvement in methane storage capacity would result. 

A comparison of the methane coverage for pure carbons having high inter- 
nal porosity with substantial pore volume contributed by pores of less than 
200 A radius (C3, C6, and C8) with those carbons having virtually no pores 
smaller than 200A radius (CIO, C15, and C16) indicate no significant differ- 
ence. Thus, micropore structure appears to contribute no additional benefits 
to methane adsorption. It should be noted that Carbons C5, C9, and Cll 
through C14 cannot be included in this evaluation because of their substantial 
inorganic content, see Table 2. 

A comparison of methane coverage versus pH for the carbon samples indi- 
cates a slight positive relationship with the most basic carbons having a high 
methane coverage. Figure 5 plots methane coverage from Table 15 versus pH 
values from Table 2. However, here too, no significance can be assigned to 
this observation since the pH contribution of the inorganic materials present 
in each sample was not determined. 

E. Conclusions 

The following conclusions can be drawn from the Task 1.1 results. 

1. Activated carbons proved superior to zeolites in this program, both on a 
grams/gram and grams/liter basis. 

2. The best performance, defined here as methane delivered per liter of 
storage, was turned in by North American Carbon G210 crushed to a packing 
density of 0.58 g/cm . 

3. High surface area and high packing density are desirable attributes for 
candidate sorbent materials. 

4. Results indicate that there is room for improvement in the performance of 
even the best carbon evaluated. 
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Figure 5. RELATIONSHIP BETWEEN COVERAGE AND pH 
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There is hope that an adsorption storage system operating at a maximum 
pressure of 3.6 MPa can be developed having a storage capacity equal to 
compression storage at 14 MPa (2000 psig). 


Table 15. METHANE COVERAGE FOR CARBON SAMPLES AT 
3.6 MPa OF METHANE 


Sample 

Designation 


Coverage in 
Molecules/cnr x IQ 1 


C3 

2.59 

C4 

2.96 

C5 

2.85 

C6 

1.99 

C7 

2.72 

C8 

3.53 

C9 

2.78 

CIO 

3.12 

Cll 

2.82 

C12 

2.70 

C13 

2.03 

C14 

1.78 

C15 

2.62 

C16 

2.83 

C17 

3.38 


Task 1.2 Modeling Effort 

The purpose of this modeling effort is to develop a simplified model 
suitable for predicting weight, capacity, and vehicle range as well as esti- 
mate the change in each parameter required to obtain a methane storage cap- 
ability for adsorption systems consistent with the range, weight, and fuel 
tank capacity of conventional (gasoline fueled) and high-pressure CNG 
vehicles. Although a rigorous analytical model was beyond the scope of this 
work, the following results are sufficient to provide a semi-quantitative 
comparison of compression storage and adsorption storage. 

A. Calculations to Determine the Weight of a Storage Tank 

Considerable developmental work has been performed to reduce the weight 
of D.O.T. approved high-pressure gas cylinders for vehicular applications. 
Aluminum-fiberglass composite cylinders having an internal volume of 38 £, an 
operating pressure of 21 MPa (3000 psig), a burst pressure of 52 MPa, and a 
mass of only 28 Kg (62 lbs) are currently available from commercial suppliers. 
However, in order to accurately assess the performance of a low pressure 
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adsorption storage system, it would be necessary to determine the weight of a 
storage tank optimized for operation at a maximum pressure of 3*6 MPa (500 
psig). Unfortunately, at the present time, there is little interest in weight 
reduction programs for cylinders operating in this pressure range. As a con- 
sequence, the 28 Kg cylinder was used as a storage container for both the high 
pressure model and the adsorption model systems, 

B, Model Automobile 

A state-of-the-art compact class automobile was chosen having the speci- 
fications listed in Table 16, These are the specifications for a gasoline 
powered Ford Tempo GL 5-speed, This size automobile was chosen for a model 
since its performance characteristics are expected to be more sensitive to 
changes in weight and fuel composition than heavier vehicles. 


Table 16. MQDEL AUTOMOBILE 


Inertial Weight 
Frontal Area 
Drag Coefficient 
Road Horsepower at 50 mph 
Aerodynamic Drag at 50 mph 
Engine Displacement 
Compression Ratio 
Power (SAE Net) 

Trunk Space 
EPA Highway 
EPA City 


2750 lbs (1250 Kg) 

20.6 ft 2 (1.91 M 2 ) 

0.36 

11.5 hp 

6.0 hp 

140 cu. in. (2301 cc) 
9.0:1 

89 bhp at 4700 rpm 
13 ft 3 (0.37 M 3 ) 

41 mpg 
27 mpg 


Road horsepower requirements and miles per gallon on gasoline were calcu- 
lated for a constant speed of 50 mph (80 km/hr) as a function of vehicle 
weight using the following assumptions: 

1 • The EPA highway mpg approximates the fuel economy at a constant 50 mph 
(80 km/ hr) 

2. Aerodynamic drag is unchanged by changes in vehicle weight 

3. Fuel consumption is related to road horsepower and engine efficiency 
requirements. 

The power required to maintain a constant velocity (P w ) is the sum of the 
power required to overcome rolling resistance (P^q) including drive train 
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losses and accessory requirements plus the power required to overcome aero- 
dynamic drag (P^)« 


P w ■ P R0 + P L 


( 8 ) 


P R0> in turn, is proportional to the product of the velocity (V), the coeffic- 
ient of rolling resistance (C^), and the vehicle weight (W). 


Pro = c R * w * v 


(9) 


Reference 9 supplies all of the necessary information required to calcu- 
late changes in road horsepower requirements as a function of changes in 
weight for our model automobile. Table 17 lists the results of these calcula- 
tions. 


Table 17. ROAD HORSEPOWER REQUIREMENTS AND MPG AS A FUNCTION OF 
VEHICLE WEIGHT AT A CONSTANT 50 MPH 


Vehicle 


Weight, lbs 

Road HP 

MPG 

2000 

10.4 

47 

2500 

11.5 

41 

3000 

12.6 

37 

3500 

13.7 

34 


The expected EPA city mpg was also calculated as a function of vehicle 
weight. The approach taken here is similar to that used by Kukkonen in 
Reference 10. Table 18 lists the results. 


Table 18. EPA CITY MPG AS A FUNCTION OF VEHICLE WEIGHT 


Vehicle 

Weight, lbs MPG 


2000 33.7 

2500 27 

3000 22.5 

3500 19.3 


C. Fuel Consumption for Dual-Fuel Mode 

Before the operational range of our model vehicle can be calculated, it 
is necessary to determine the fuel consumption expected during natural gas 
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operation# Since a gasoline fueled spark ignition (SI) engine can be made to 
run on natural gas with only minor modifications, it has become standard 
practice in the U.S.A. and elsewhere to perform dual-fuel conversions. Such 
conversions allow the driver to select the fuel desired and to change from one 
fuel to the other at will. Such conversions do not take advantage of all of 
the positive attributes of natural gas as an SI engine fuel but, since such 
conversions are common, the decision was made to include a dual-fuel mode in 
our model. 

Work carried out by the University of British Columbia allows some quan- 
tification of the fuel consumption for gasoline engines converted to natural 
gas fuels. ^ This work indicates that the overall operating efficiency of a 
vehicle equipped with a natural gas-fueled converted gasoline engine is, at 
best, 12% better on natural gas than the same engine on gasoline, provided 
modifications are made to the timing curve to advance the spark when operating 
on natural gas. A more detailed discussion of engine efficiency in the dual- 
fuel mode can be found in Part 1 of Appendix C. 

D. Dedicated Engine Mode 

Natural gas exhibits many excellent qualities as an internal combustion 
(IC) spark-ignited engine fuel. Chief among these qualities is the high 
octane number of 130 RON. However, to take advantage of this high octane 
characteristic, it is necessary to use high compression ratios which are not 
compatible with gasoline operation. 

An engine designed to make optimum use of natural gas as a fuel is 
classed as a dedicated engine. Such an engine could be expected to be 25% 
more efficient on natural gas than the same displacement size engine running 
on gasoline. A more detailed discussion of the efficiency of a dedicated 
engine can be found in Part 2 of Appendix C. 

E. Energy Density of Pressurized Storage 

In addition to knowing the efficiency and, therefore, the range per unit 
mass of fuel, it is also necessary to know the total amount of fuel stored on 
board. Because the term "pipeline quality" natural gas does not carry with it 
a definitive gas composition, it is necessary, at this point, to define terms. 
For purposes of this model, natural gas will have an average molecular weight 
of 16.0 and a lower heating value of 11,800 Kcal/kg (21,250 Btu/lb). By 
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placing the average molecular weight equal to that of methane, we can apply, 
without further adjustment, the data gathered for methane adsorption on car- 
bons. This leads, however, to underestimating vehicle range by a factor equal 
to the ratio of the actual average molecular weight of a specific gas sample 
to that of methane. 

The amounts of natural gas (methane) that can be stored per liter of 
volume at various pressures is summarized in Table 19. 


Table 19. RELATIONSHIP BETWEEN PRESSURE AND METHANE STORAGE 
CAPACITY AT 23 °C FOR AN EMPTY CYLINDER 



Pressure 


Energy Density, 

P si S 

atmospheres 

M Pa 

Grams/Liter 

500 

35 

3.6 

24 

1000 

69 

7.0 

50 

1500 

103 

10.4 

79 

2000 

137 

13.9 

108 

2500 

171 

17.3 

135 

3000 

205 

20.8 

160 

We can now utilize 

this table 

to determine the 

fuel on-board given the pres- 

sure and volume of 

the containment vessels. There is one additional correc* 

tion that must be 

applied to 

an actual system. 

State-of-the-art dual-fuel 


systems cycle between a maximum tank pressure of 16.7 to 20.8 MPa and a 
minimum tank pressure of 0.3 to 0.4 MPa depending on manufacturer. Thus, of 
the 108 grams/liter stored at 13.9 MPa, only 106 g/£ is delivered by cycling 
the system from 13.9 MPa to 0.3 MPa. However, for purposes of calculation, 
the higher value of 108 g / i will be used to offset somewhat the error induced 
by using 16 as the average molecular weight of natural gas. 

F. On-Board Storage 

The calculation of the "fuel on-board" for an adsorption system is more 
complex than for an empty cylinder. Gas can be stored both by adsorption on 
the substrate and by compression in void spaces. Furthermore, the quantity of 
methane adsorbed is not a linear function of pressure but follows a Langmuir- 
type isotherm. Rather than address these problems from a theoretical point of 
view, actual experimental values for the best carbon identified in Task 1.1 
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were used. North American Carbon G-210 crushed to a packing density of 0.58 
g/cc was loaded into a cylinder and cycled between 3.6 MPa and 170 kPa (500 
psig and 10 psig respectively). About 58 grams of methane were delivered per 
liter of carbon during the discharge cycle. It is premature to conclude that 
the 58 g / i of methane delivered by the best carbon represents a maximum value 
for absorbent systems in general. 

G. Additional Weight of Gas Storage Hardware 

In order to determine the performance and range on natural gas of a 
vehicle operating in the dual-fuel mode, it is necessary to know the addi- 
tional weight added to the vehicle by the gas storage system. Table 20 
summarizes the additional weight of one, two, and three tank storage systems 
based on 38-liter, internal volume, aluminum-fiberglass composite cylinders. 
The 28-Kg high pressure cylinder previously described is used to calculate the 
weight of the various storage systems. The packing densities of the carbon 
adsorbents are set at 0.58 g/cc, which represents the most dense carbon 
evaluated in this program. A single 38)1 tank with manifold and bracketing 
hardware would occupy about one-third of the available trunk space in our 
model vehicle. 


Table 20. APPROXIMATE MASS OF STORAGE SYSTEMS IN Kg 


Pressurized Storage 

1 Tank 

2 Tanks 

3 Tanks 










Fuel 

6 

12 

18 

Tanks 

28 

56 

84 

Brackets 

10 

20 

30 

Conversion Kits 

25 

25 

25 

Total 

69 

113 

157 

AdsorptiQn Storage 

1 Tank 

2 Tanks 

3 Tanks 

Fuel 

3 

5 

7 

Tanks 

28 

56 

84 

Brackets 

15 

30 

45 

Manifold 

30 

30 

30 

Carbon 

22 

44 

66 

Total 

98 

165 

232 
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H. Vehicle Range and Fuel Efficiency 

Having determined the weight and storage capacity of the various storage 
systems and the relative fuel efficiencies of the SI engine operating on gaso- 
line and on natural gas in both dual-fuel and dedicated modes, it is now pos- 
sible to combine these relationships with the relationship between total 
vehicle weight and vehicle fuel economy to calculate both gasoline equivalent 
mpg and range on natural gas for the various storage options. Table 21 com- 
pares fuel efficiency on gasoline and the estimated range on natural gas for 
the various storage options for a model vehicle operating in the dual-fuel 
mode on the EPA-city cycle. Table 22 is a similar comparison for the constant 
speed approach. Results in Tables 21 and 22 assume that pressurized storage 
system cycles between 20.8 MPa (3000 psig) and 0.3 MPa (25 psig) while adsorp- 
tion storage systems cycle between 3.6 MPa (500 psig) and 0.2 MPa (10 psig). 
Range on natural gas is based on the results of the study discussed in 
Appendix C. 

The projected range of the model vehicle equipped with a dedicated engine 
would be about 10% longer than those listed for various storage options in 
Tables 21 and 22. 

I. Discussion of Results 

Fuel Economy : The results of our model indicate that a significant 

x penalty is paid, up to 15%, in the form of a loss in fuel economy on the EPA 
city cycle when a dual-fuel capability is added to the base automobile. On 
the other hand, there is little additional performance loss with adsorption 
storage over that experienced with pressurized storage. This is an important 
result since one of the past objections to adsorption storage has been the 
added weight of the sorbent bed. This has added significance since our model 
calls for heavier brackets and manifolding hardware with adsorption storage 
than with pressurized storage. In other words, there is no significant 
difference between pressurized storage and adsorption storage in the fuel 
economy loss experienced in the dual-fuel mode. 

Range : The range of a vehicle is directly related to the fuel stored on- 

board and it is generally accepted in alternatively fueled vehicle circles 
that for an alternatively fueled vehicle, regardless of the propulsion system, 
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Table 21. ESTIMATED FUEL ECONOMY AND RANGE ON THE EPA CITY CYCLE 
FOR VEHICLE OPERATING IN A DUAL-FUEL MODE 


Gasoline Fuel 

Storage System Economy* 

Gasoline Only 11.5 (27) 

Pressurized Storage (20.8 MPa) 


Range on 

Natural Gas Only** 


1 Tank 

10.9 

(25.4) 

95 

( 59) 

2 Tanks 

10.5 

(24.6) 

182 

(113) 

3 Tanks 

10.1 

(23.7) 

262 

(163) 


Adsorption Storage (3.6 MPa) 






1 Tank 

10.7 

(25.1) 

32 

( 

20) 

2 Tanks 

10.2 

(23.9) 

63 

( 

39) 

3 Tanks 

9.7 

(22.8) 

93 

( 

58) 


* Kilometers per liter (miles per gallon). 

** Kilometers (miles) . 


42 



Table 22. ESTIMATED FUEL ECONOMY AND RANGE AT A CONSTANT SPEED 
OF 80 Km/hr IN THE DUAL FUEL MODE 

Gasoline Fuel Range on 


Storage System 


Gasoline Only g 18 

Pressurized Storage (20.8 MPa) 

1 Tank 17 

2 Tanks 17 

3 Tanks 16 

Adsorption Storage (3.6 MPa) 

1 Tank 17 

2 Tanks 16 

3 Tanks 16 


* Kilometers per liter (miles per 

** Kilometers (miles) . 


Economy* Natural Gas** 


(41) 



(40) 

148 

( 92) 

(39) 

290 

(181) 

(38) 

428 

(266) 


(39) 

51 

( 

32) 

(38) 

101 

( 

63) 

(37) 

150 

( 

93) 


gallon) . 
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to make significant market penetration it must have a range of at least 100 
miles between refuelings on the EPA city cycle. Although this is less criti- 
cal for a dual-fuel vehicle, it is absolutely necessary for a dedicated 
vehicle. 

As can be seen in Table 21, our model predicts that this can be achieved 
with a pressurised storage system operating at 20.8 MPa (3000 psig) equipped 
with two 38 % storage tanks. On the other hand, the state-of-the-art adsorp- 
tion system with three 38£ tanks can only manage 58 miles (93 km). It is 
possible that future research and development can boost the storage capacity 
of sorbent beds based on carbon to the point where a 100 mile range on the EPA 
city cycle is feasible with three tanks. Also, it should be recognized that 
this model assumes that storage is limited to the space available in the auto- 
mobile trunk and assumes cylindrical shaped containers, the maximum possible 
being three 38fc tanks. If, however, storage volumes in excess of 180& (about 
five 38 tanks) can be made available either by locating tanks in other unused 
spaces, by redesigning the automobile, or by using more efficient geometrical 
shapes, a state-of-the-art adsorption system could achieve a range in excess 
of 100 miles on the EPA city cycle at the expense of an additional weight 
increase without further improvements in the sorbent beds. It is unlikely, 
however, that an adsorption storage system will ever be designed with the 
energy density of a pressurized storage system charged to 20.8 MPa. 
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TASK 2. LITERATURE SURVEY AND ADVANCED STORAGE MEDIUM EVALUATION 


The objective of this task was to survey the available literature and 
identify materials and concepts with the potential to store methane at low 
pressures* The survey concentrated on the concepts of clathration, encapsu- 
lation, and dissolution* A 33-page evaluation of clathrate storage concepts, 
complete with its own bibliography, is included as Appendix D. A 27-page 
discussion of methane storage by dissolution is included as Appendix E* A 
summary of the findings is provided below* 

A. Clathration of Methane 

The terra clathrate was originally defined as a compound " in which two or 

more molecular components are associated without ordinary chemical union but 

through complete enclosure of one set of molecules in a suitable structure 
12 

formed by another* Thus, the term clathrate implies guest molecules caged 
within spaces formed by the crystal structure of the host. The phenomenon is 
dependent primarily upon the relative size of the holes in the host crystal 
structure and size of the guest molecules. 

Although there is considerable work in progress toward developing the 
know-how for tailoring new clathrate hosts for specific guest molecules, most 
of the work is oriented toward biological and catalytic applications such as 
enzymes and ionic species* Practically no references have been found for 
tailoring the clathrate host molecules for small and relatively inert mole- 
cules such as methane or the inert gases; thus, tailor-making host molecules 
for clathrating methane still requires a largely empirical approach in devel- 
opment. Nonetheless, there are a few precepts available to guide such a 
development: 

• Clathrate formulation depends on the ability of the host molecules to 
form a crystal habit which has cavities large enough to accommodate the 
specified guest molecules* 

• The crystal habit formed in the presence of the guest molecule is usually 
different from the normal crystal structure found in the absence of the 
guest* 

• The primary prerequisite for the host is hydrogen-bonding or other com- 
plex forming ability while a secondary prerequisite is a molecular geom- 
etry which will crystallize in the presence of the guest to a structure 
with cavities sized to accommodate the guest* 
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Stability is of major importance since it controls not only whether the 
methane can be contained at reasonable pressures under ambient conditions but 
also the conditions under which the methane can be released. The hydrate of 
methane, for example, is not stable enough since it requires a pressure greater 
than 28 MPa to exist under ambient temperature conditions. On the other hand, 
the hydroquinone/methane clathrate is too stable. Thus, although it can be 
stored at ambient temperature and pressure, it requires an elevated temperature 
or the introduction of a solvent to release the methane. Obviously, a com- 
promise is in order so that the clathrate will be stable at ambient tempera- 
ture, at a reasonable pressure of 1.5 to 3.6 MPa, so that methane release can 
be achieved and controlled by pressure reduction. 

The parameters controlling stability include: 

• Hydrogen bonding power of the groups through which the crystallization 
occurs. 

• Geometry and symmetry of the host molecule as it affects the structure of 
the crystal habit formed. 

The capacity of a clathrate for methane, on the other hand, depends on 
two factors: 

• How many host molecules are required to provide one guest "cage” (i.e., 
the unit cell) 

• Molecular weight of the host molecules. 

The structure of the unit cell and, therefore, the number of host molecules 
per guest molecule is dependent in a very complex way on various crystallo- 
graphic factors and is beyond prediction at this time except perhaps for cer- 
tain "hexamer" clathrates which tend to give a 3:1 host to guest ratio. Thus, 
molecular weight may be the factor of greatest importance in determining 
capacity and is suitable for preliminary evaluation. 

It may be possible that some higher molecular weight molecules with com- 
plex symmetries form clathrate structures in which the host/guest ratio is 
significantly less than the 3:1 found for the hexaraer types. Others may form 
structures with cavities sufficiently large to accommodate two or more methane 
molecules as in the case of the urea tunnel adducts. In the latter case, dif- 
fusion would be the primary barrier to decomposition and the system would 
function in a manner analogous to that of absorption in a zeolite. Most of 


46 



the clathration entities forming larger cavities (and there are a lot of them 
including Dianin's compound, the Werner complexes, deoxycholic acid, urea and 
thiourea cyclodextrin, etc.) have been evaluated only with larger guest mole- 
cules (usually the solvent used) at atmospheric pressure and not with perma- 
nent gases such as methane under elevated pressures. For example, urea adduc- 
tion of the homologous n-paraffin series has been extended downward to include 
propane and butane which are stable at atmospheric pressure only at subzero 
temperatures. However, it is possible that adduction could be extended to 
include methane and ethane as well at ambient temperatures and elevated pres- 
sures. 

It is apparent from the above discussion that any future work toward the 
development of clathrate systems for on-board storage of natural gas will be 
highly empirical in approach, although guided by the precepts outlined above. 
Two as yet unreported possibilities are a urea-methane adduct and an 
acetamide-methane clathrate. 

1 . Urea/Methane Adduct 

The methane/urea adduct, if it exists, probably will show a mole ratio of 
about 1:3 or a weight ratio of about 16:180 » 0.089 g CH^/g of urea. This is 
of the same magnitude as that for activated carbons but is not as good as the 
methane hydrate (0.155 g/g H 2 0). 

2. Acetamide/Methane Clathrate 

o 

Acetamide CH 3 CCNH 2 MW = 59 MP - 82°C 

Methane clathrates of acetamide have never been demonstrated but are good 
candidates because of good hydrogen-bonding ability of the amide group and low 
molecular weight. If a host/guest ratio of 3:1 is achieved, the weight ratio 
would be 0.09 g CH^ per g acetamide. The possible methods of preparation in- 
clude crystallization from a solvent or condensation from a vapor. 

Conclusions: 

As reported above, the highest mass ratio of guest to host is found in 
methane hydrates, 0.155 g/g. However, methane hydrate is not stable under 
ambient temperature conditions at pressures below 28 MPa. The best mass ratio 
that can be expected from other proposed but as yet unidentified methane clath- 
rates is about 0.09 g/g. This is slightly less than the mass ratio already 
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available with state of the art carbon based adsorption systems. Also, since 
a clathrate-based storage system is likely to be relatively more difficult to 
cycle than an adsorption-based storage system, further work with clathration- 
based systems for on-board storage is not recommended. However, clathration- 
based systems may ultimately prove useful for stationary applications where 
weight and volume are not critical. Therefore, recommendations for future 
research with clathrates is presented at the conclusion of Appendix D, despite 
its lack of utility as an on-board storage concept. 

B. Dissolution of Methane 

In general, the data conforms well to the classical views of Hildebrand J 
on regular solutions, i.e., solutions M in which orienting and chemical effects 
are absent and in which the distribution and orientations are random..." In 
other words, methane behaves in most solvents, including the rare gases, as 
non-polar molecules subject primarily to van der Waals dispersion forces. The 
primary parameters controlling solubility in regular solution theory are the 
molal volume and the Hildebrand solubility parameter relative to that of 
methane. The solubility parameter or cohesive energy density is defined as: 


S 


( 



1/2 


( 10 ) 


where: 

Afl? = Heat of vaporization at solution temperature 

Vj = Molal volume of Component 1. 

Sipce the solubility parameter of methane is at the lower end of the 
scale (S - <6), the best solvents for methane are the perf luorocarbons, the 
silicones, and the lower aliphatic hydrocarbons such as propane. The best 
solvent for methane so far found is propane, in which methane is soluble to 
the extent of 0.15 mole fraction or 0.063 lb CH^/lb solvent. However, this is 
a lower capacity on a lb/lb basis than can be expected from adsorption on 
activated carbon at this pressure. 

On a mole fraction basis, the solubility of methane in octamethyl cyclo- 
tetrasiloxane (0.32) and in perfluoro p-heptane (0.28) is greater than the 
propane but on a gravimetric basis the solubilities are much lower (0.025 and 
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0.016, respectively) due to the high molecular weights of these solvents* 
Nevertheless, it is possible that lower molecular weight per fluorocarbons and 
silicones could exhibit acceptable methane solubilities on a mass ratio basis. 
However, information regarding methane solubility in these solvents is not 
available in the literature. 

In the aliphatic hydrocarbon series, the solubility of methane decreases 
with increasing molecular weight as the result of the solubility parameter in- 
creasing with increasing molecular weight. This trend of decreasing methane 
solubility with increasing molecular weight continues up to about (n- 
hexane). Above Cg, however, the solubility of methane on a mole fraction 
basis increases with increasing molecular weight to 0.32 for C 30 (squalane) . 
Unfortunately, the solvent molecular weight increases more rapidly than solu- 
bility so that the solubility on a gravimetric basis continues to decrease 
with increasing molecular weight. These results are based primarily upon data 
for straight chain hydrocarbons, squalane being the exception with 6 methyl 
groups oh a straight chain of 24 carbons. However, the solubility of methane 
in isobutane (2-methyl propane), 0.061 g/g, and in neo pentane (2,2-diraethyl 
propane), 0.057 g/g, suggest that the decrease in methane solubility on a 
gravimetric basis with increasing molecular weight may be less severe for 
highly branched paraffin hydrocarbons. Unfortunately, data on the solubility 
of methane in highly branched high molecular weight paraffins is not available. 

The solubilities of methane in aliphatic alcohols appear to be greater, 
at the same solubility parameter value, than for other solvents having no 
hydroxyl groups. This suggests that hydrogen bonding may contribute in some 
as yet unidentified manner to enhancing methane solubility. This further 
suggests that it may be possible to enhance the adsorptive capacity of solid 
adsorbents by providing a multiplicity of hydroxyl groups on the surface. 

The solubility data for the lower paraffin hydrocarbons (e.g., propane) 
actually represent the concentration of methane in the liquid phase of a mix- 
ture in equilibrium with a vapor phase containing both methane and solvent 
vapor. If the equilibrium concentration of methane in propane at 25°C and 3.6 
MPa total pressure is corrected to 3.5 MPa partial pressure of methane, the 
solubility would increase from 0.063 g/g to 0.11 g/g but the total system 
pressure would also increase. Likewise, the solubility of methane in liquid 
ethane at 25°C and a methane partial pressure of 3.6 MPa is 0.35 g/g, however, 
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system pressures exceeding 7 MPa would be required. If we could tie down 
ethane or propane as ethyl or propyl groups and still maintain their solvent 
capacity for methane, gravimetric capacities similar to those of the best 
adsorbent systems might be achieved. 

Conclusions 

No high molecular weight solvent has been identified which can dissolve 
enough methane to be competitive with adsorption storage. On the other hand, 
no method has been identified to prevent low molecular weight highly volatile 
solvents such as propane, in which methane is highly soluble, from vaporizing 
along with the methane during system discharge. Although this is not consid- 
ered a serious problem with propane, it does present a significant environmen- 
tal hazard when low molecular weight fluorocarbons and silicones are considered. 
The possibility of tying low molecular weight solvents to polymer backbones, 
however, is at least within the realm of speculation, and warrants some future 
consideration. 
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TASK 3. FUTURE RESEARCH AND DEVELOPMENT RECOMMENDATIONS 


Theoretical and experimental results obtained in this program clearly 
indicate the potential of low-pressure natural gas storage systems (based on 
carbons as the storage medium) for vehicular applications. As the results 
show, by adsorption on carbons commercially available today, it is possible to 
store at 500 psi approximately 70% of the volume of methane that can be stored 
at 2000 psi in a compressed (no adsorbent) natural gas system of the same 
volume, and about 45% of that stored at 3000 psi. 

The cost of compressors capable of compressing natural gas to pressures 
greater than 2500 psi contributes a significant amount to the total cost of 
natural gas-fueled vehicle fleets. Consequently, lowering the compression 
requirements could significantly impact the total system cost and make natural 
gas a more attractive transportation fuel. 

The potential of low-pressure storage systems has been demonstrated, al- 
though the range of vehicles operating on adsorption systems is less than for 
high pressure systems. Consequently, there is need for additional research, 
development, and optimization before these low-pressure systems can be imple- 
mented. Following is a list outlining the areas that must be resolved. 

1) There is a need to continue research with carbons to improve their stor- 
age capacity. The likely approach to be taken is the manipulation of the 
surface chemistry of the carbon. As was reported in Task I, certain of 
the carbon samples gave a higher degree of coverage of methane than 
others, most notably those with high pH values on the acid-base test. An 
understanding of the specific surface groups which enhance methane 
adsorption can lead to a method of increasing methane storage capacity by 
increasing the population of the desirable surface group. 

2) During the adsorption of methane on carbon, a significant amount of heat 
is generated. Although some work has been performed with systems using 
5 liter storage cylinders, no one has addressed the problem of heat dis- 
sipation in larger cylinders. 

3) Golovoy^ has observed that hydrocarbons larger than are strongly 
adsorbed on high surface area carbons during the charging cycle to the 
point that they do not desorb on discharge. It is necessary that the 
impact of these materials upon cycling efficiency be determined and 
several methods for removing them from the fuel prior to loading be 
evaluated. 

4) It is expected that the use of low pressure adsorption systems in place 
of high pressure storage systems will reduce the capital investment re- 
quired for a refueling station but may increase the cost of the vehicular 
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storage system. It is necessary to conduct an economic evaluation to 
quantify the expected savings which will be realized by reduced com- 
pressor costs so that an upper limit can be placed on the cost of the 
vehicular adsorption storage system. 

5) To date, most adsorption system experiments have been conducted at about 
25°C. However, real world methane storage systems may be exposed to 
environmental extremes of from -30° to -H>0°C with daily temperatures 
swings of as much as ±20°. It is necessary to determine the impact upon 
both performance and safety of daily and annual temperature variations. 

6) Vibrational settling of adsorbent may present a problem in actual vehic- 
ular systems. This problem, if it exists, needs to be identified early 
enough for corrective measures to be taken prior to actual full-scale 
demonstration tests. 

The above issues need to be resolved in order that enlightened full-scale 

demonstration tests can be undertaken. 
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PORE SIZE DISTRIBUTION FOR CARBON SAMPLES 
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APPENDIX C 


SUPPORTING CALCULATIONS FOR TASK 1.2 WORK 

I • Estimation of the Fuel Economy of a Gasoline Engine Converted to Run on 
Natural Gas In a Dual Fuel Mode 


• Model vehicle achieves 17 km/jt (40 mpg) or 22.9 km/kg (6.54 miles/lb) on 
gasoline 

• Lower heating value of gasoline is 10.5 kcal/gram (18,900 Btu/lb)* 

• Lower heating value of natural gas is 11.56 kcal/gram (20,800 Btu/lb)* 

Work carried out by the University of British Columbia (UBC) has provided 
some quantification of the fuel economies for vehicles operating in a dual 
fuel mode on gasoline or natural gas. The UBC study documented an average 
12.5% increase in operating efficiencies for vehicles operating on natural gas 
over the same vehicle operating on gasoline on the basis of BTU input. This 
translates into a 23% increase in fuel economy on a mass to mass basis. In 
other words, the UBC results suggest that our model vehicle which achieves 
22.9 km per Kg of gasoline will travel 28 km per Kg of natural gas. 

The proposed explanations for such an improvement in fuel economy were a 
higher thermal efficiency for the engine operating cycle due to more complete 
combustion of the natural gas-air mixture and a leaner stoichiometric ratio in 
the natural gas mode. However, the performance of the older carbureted engines 
used in the UBC study are not representative of the high efficiency, electronic 
fuel injected engine used in the model vehicle chosen for this study. Conse- 
quently, most of the 12% efficiency improvements, on a Btu basis, would not be 
achieved with a state-of-the-art engine. Not because of a decrease in perfor- 
mance on natural gas, but rather a result of improved performance and fuel 
economy on gasoline. 

For the purposes of the modeling work in Task 1.2, it was assumed that 
1 kilogram of natural gas delivered about 12% greater range than 1 kilogram of 
gasoline, 10% as a consequence of the higher energy content of natural gas and 
2% as the result of improved engine efficiency. 

To conduct a more meaningful analysis of the potential efficiency improve- 
ment or vehicle range (beyond that performed here) would require the following: 

m Modeling of natural gas-air and gasoline-air cycles for selected design 
specification of IC engines and numerical simulation of the performance 
within a defined range of operating conditions 
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• Selection, analytical modeling, and simulation of a reference driving 
cycle in terms of determined natural gas-air and gasoline-air cycle 
performance 

0 Comparison analysis of driving cycle performance for both natural gas and 
gas fuels carried out in terras of BSFC or distance unit of fuel mass. 

II Estimation of the Fuel Economy of a Dedicated Natural Gas Engine 

In general, natural gas exhibits excellent characteristics for use as an 
internal combustion spark-ignited engine fuel. The clearly identifiable advan- 
tages and disadvantages are as follows: 

Advantages 

1. High caloric value (lower heating value is 10% higher than gasoline) 

2. High octane number (130 RON) 

3. Good ignitability of natural gas-air mixture allows low equivalence 
ratios 3 ^ 

4. Simple oxidation reactions leading to complete combustion 

5. Effective carburetion due to gaseous phase. 

There are also certain negative natural gas characteristics: 

Disadvantages 

1. Relatively low laminar flame speed within quiescent combustion chamber 

2. Lack of process of vaporization both during mixing and admission elimin- 
ates the cooling effect and accordingly decreases cylinder volumetric 
efficiency. 

Advantages 1, 4, and 5 can be made use of in a careful conversion of a 
gasoline engine. The first disadvantage can also be partially corrected by 
advancing the ignition and also to varying degrees by the type of engine that 
is selected for conversion. 

More complete utilization of advantageous characteristics of natural gas 
as an IC engine fuel as well as effective alleviation of the negative charac- 
teristics will require a dedicated natural gas-engine design. The most evi- 
dent modification should Involve: 

0 Increase of compression ratio (up to 14 to 1) in order to maximize 
thermal efficiency of cycle by extensive utilization of high octane 
number of natural gas 


C-2 



• Application of "fast burn" combustion system in order to alleviate 
disadvantage of slow flame propagation of natural gas-air mixture and 
enhance cycle thermal efficiency by bringing cycle conditions closer to 
those characterizing a const ant- volume cycle 

• Application of some sort of supercharging to alleviate the deterioration 
of volumetric efficiency caused by the lack of cooling upon carburetion. 

Figure C-l defines the functional relation between theoretical thermal 
efficiency of an internal combustion-spark ignition engine and design compres- 
sion ratio developed for gasoline-air mixtures with different equivalency 
ratios F r .^ The process of transition from a converted natural gas engine to 
a "dedicated”* natural gas engine can be represented in this figure by the 
transition from Point A to Point B. Point A represents an engine with a 
compression ratio Tq » 9 typical for a modern gasoline engine operating on 
slightly over-stoichiometric mixture with F^ = 1.1 which seems to be approp- 
riate for assuring the satisfactory combustion process in a gasoline engine 
converted to natural gas fuel. Although lean methane-air mixtures have good 
ignitability , their poor flame propagation characteristics seem to prevent 
effective burning in a conventional gasoline mixture combustion chamber. 
Satisfactory fast combustion of natural gas mixtures in such conditions will 
call for an ignition advance and for a rather enriched equivalency ratio. A 
dedicated design can be provided with a compression ratio r^g = 11.5 and a 
combustion system intentionally promoting faster combustion, therefore capable 
of accepting leaner mixtures. Point B can be then quite realistically extrap- 
olated as to be located on curve F^ = 0.8 and abscissa r^Q " 11.5. This ex- 
trapolation suggests a possible corresponding thermal efficiency gain An' « 8%. 
An' can still be increased by supercharging allowing for approximately a 5% 
additional increase in theoretical thermal efficiency of an IC engine. The 
total efficiency increase for a dedicated natural gas IC engine will then be 
approximately 13% over a natural gas converted version, or if the gasoline 
engine is taken as the baseline, the potential efficiency improvement of a 
natural gas dedicated design would be about 25%. 

To a first approximation, the ranges of the model vehicle described in 
Task 1.2 and listed in Tables 21 and 22 for natural gas operation can be 
extended by 13% if a dedicated engine is substituted for a converted gasoline 
engine. 
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APPENDIX D 


EVALUATION OF CLATHRATION COMPOUNDS AS A MEANS OF 
STORING NATURAL GAS 

CLATHRATES 1-16 


1 • Background 

Interest in inclusion compounds in recent years has been found in two 
major areas. In the firht and most recent, the binding or complexing of guest 
species by unimolecular hosts, in solution , has received much attention, par- 
ticularly in the fields of biology and catalysis. The second relates to the 
study of crystalline inclusion compounds or clathrates which may be subclas- 
sified as: 

• The cage clathrates in which the guest molecules are imprisoned in dis- 
crete closed cavities or cages in the host crystal and include the hydro- 
quinone and ice clathrates 

• The channel type in which the guest species are accommodated in continu- 
ous channels in the crystal such as the urea and thiourea adducts 

• The layer type* in which the guest species are situated between layers 
(e.g., with graphite and certain clays). 

In addition, molecular sieves and, in particular, the zeolites are sometimes 

classified "as clathrate or inclusion hosts since they possess discrete cages 

and channels, and inclusion depends on the molecular dimensions of the guest 

17 18 

molecules relative to that of the channels and cages of the host. * From 
the point of view of methane storage, the clathrate type of inclusion is of 
greatest interest. However, there may be a large overlap in concept and prin- 
ciple between inclusion in the solutions and in clathrates. Furthermore, 
inclusion in solutions may provide a means of significantly enhancing methane 
storage. 

The term clathrate was originally defined, on the basis of x-ray diffrac- 
tion structural analysis of a number of such entities by Powell and associates, 
as a compound "in which two or more molecular components are associated with- 
out ordinary chemical union but through complete enclosure of one set of 
molecules in a suitable structure formed by another."^ Two aspects of the 
clathrates are noteworthy: 

• The crystalline lattice structure of the host is usually not in its nor- 
mal crystalline mode (a”form) but in a form (usually referred to as the 
g-form) which is usually less stable in the absence of guest molecules 
than the orform. 


D— 1 



• The composition over which the clathrate is stable with respect to decom- 
position into guest and normal host can vary over a large range. 

Thus, clathrates are not stoichiometric compounds but ones in which the re- 
strictions on which molecules can become guests and the minimum fraction of 
the total number of cavities which must be filled for stability are largely 
determined by geometric considerations. 

Although the literature concerning clathrate inclusion compounds extends 
back to the early 1800’s, most of the activity in the field has occurred since 
the late 1940’s and early 1950 ’s after a basic understanding of the phenomenon 
was developed by H. M. Powell and associates. Since then, the work in the 
field has been prodigious. Thus, Bhatnager (1970)* lists 1339 references which 
is certainly not complete since this includes only 293 references on hydrates 
compared to 1458 references cited by Davidson. Nonetheless, Bhatnager’s bib- 
liography yields an interesting breakdown of the types of clathrating agents 


of greatest interest: 

• Urea and thiourea adducts 610 refs. 

• Hydrates 293 

• Phenolics including hydroquinone, 

Dianin’s compounds, and various phenols 116 

• Cyclodextrins (starches) 63 

• Werner complexes 134 

9 Miscellaneous including phosphonitrile, 

steroids, adamantane, cholesterol, 
dinitrophenol , hexamethylisocyanide 

chloride, and cycloveratril 41 

9 Reviews 82 


It is interesting to note that many of these clathrating agents exhibit a 
strong hydrogen bonding propensity (e.g., amides, water, phenols, and cyclo- 
dextrin) which is undoubtedly involved in the binding of the crystals of the 
host molecules. In the case of the Werner Complexes and other inorganic co- 
ordination systems, the binding energy for the host crystal certainly involves 
ionic as well as ion to ligand attractions. 

Since 1967, the work in the field has begun to shift toward developing 
the know-how to synthesize new host structures, a priori, on the basis of 
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molecular structure rather than on a purely empirical basis. One area of con- 
siderable interest in the biological field is the synthesis of water soluble 
host molecules, which are organic molecules that contain cavities of dimen- 
sions capable of accommodating simple ions and molecules in which the binding 
is provided by hydrogen bonding, ion pairing, metal ion to ligand attractions, 
acid-base attractions, and van der Waals forces. Hosts may contain cavities 
which result from reorganization of the molecule during complexation with a 
guest molecule or rigid cavities which exist prior to complexation 
(cavitands) • ^ 

This work is an outgrowth of the studies of the structure and properties 

of the cyclodextrins, which are naturally occurring cyclic oligomers made up 

of 6-8 glucoside units bound head-to-tail and which enclose cavities of 6, 8, 

or 10A diameter, depending on the number of units in the oligomer. The ether 

oxygens and hydrogen atoms are oriented inwardly and hydroxyls are oriented 

outwardly. Thus, guest molecules of a range of sizes from methane to much 

larger organic species can be clathrated. Similar behavior was first achieved 

with the synthesis of cyclic polyethers containing six (-CHj-CHj-O-) units 

with inwardly oriented ether groups. Furthermore, this approach has now 

been extended to include macrocycle compounds made up of substituted aromatic 

units, for which the size and shape of the cavities as well as the nature of 

21 

the inwardly oriented binding groups has been varied. 

14 21—2*1 

Independently, MacNicol and associates > have developed a synthesis 

strategy based on mimicking the "hexa-host" behavior of several phenolic clath- 
rating agents such as Dianin's compound and hydroquinone, in which clathration 
involves the formation of a hexagonal ring of hydrogen-bonded oxygen atoms 
from six phenolic hydroxyl groups: 
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Each group of six substituted phenol molecules is associated with two cavi- 
ties , one above the hexagonal hydrogen bonded planes and one below. 

They were then able to demonstrate that the clathration ability of such 
hexagonal structures could be simulated using the permanent structure of cer- 
tain hexa-substituted benzene ring compounds* Thus, two independent approaches 
starting with empirical observations of two radically different types of clath- 
rates have converged on a cyclic model for the synthesis of new clathrate hosts 

However, the guest molecules bound by both types of clathrate hosts so far 

2 1 22 

have been relatively larger than methane. * 

Thermodynamics of Clathrate Stability^ 

The thermodynamics of clathrate stabilities has been reviewed by Child 
(1964)^ and will not be repeated here. However, the gist of his analysis is 
of interest in guiding further synthesis attempts. He points out that in many 
respects clathrates are similar to equilibrium solutions in that: 

• Solute molecules are placed in cavities within the solvent 

9 The energy of interaction between host and guest is normally small and 
the entropy term is similar in magnitude to the entropy of vaporization 
of a solute from a solution which obeys Raoult’s law 

9 The guest may be said to obey Henry’s law in the sense that guest-guest 
interactions are negligible compared to the guest-host interactions. 

On this basis. Child suggests that the entropy term contribution to sta- 
bility is relatively small and that the stability depends primarily on the 
overall heat of formation of the clathrate and, therefore, on the magnitude of 
interaction between host molecules, and between the guest and the surrounding 
molecules of the host lattice. Thus, the heats of formation (AHp) for a 
number of clathrates formed from hydroquinol and water have been shown to be 
about 0.8-2. 2 of the value of two times the heat of vaporization AH vap of the 
guest molecules, i.e., AHp = 0.8 to 2.2 (2AH vap ). 

An explanation for these high ratios was given in terras of: 

9 The removal of a guest molecule from a clathrate leaves an empty cavity 

(with little interaction between host molecules across the hole), whereas 
the removal of a molecule of guest from a liquid leaves no hole 

9 The interactions between guest molecules and the wall of the cage may be 
greater than between the two guest molecules in the liquid state 
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® The coordination number of the guest molecule may be larger in the cavity 
than in the liquid state. Thus, the vaporization of the pure guest may 
be visualized as a two-step process in which the vaporization step leaving 
empty holes in the remaining liquid has an enthalpy increase of 2 aH v , 
and the subsequent collapse of the holes an enthalpy decrease of “AH vap 
for a net change of AH . In the case of clathrates, on the other hand, 
the vaporization enthalpy change of 2 AH is counterbalanced not by the 
collapse of the holes remaining, but by the enthalpy change in converting 
the 8 structure to the a-form which is often much less than AH p . Fur- 
thermore, the heat of vaporizing the guest molecule from the clarhrate 
structure may actually be greater than the heat of vaporization of the 
pure liquid because of the interaction between guest and host by an 
amount depending on the energy of interaction (which depends on the chem- 
ical nature of the guest and host) and the coordination number of inter- 
action (which will depend on the size and shape of the cavity relative to 
that of CH^). 

Since methane is non-polar in nature, the host-guest interactions 
are likely to be limited to van der Waals forces. Thus, the most impor- 
tant host-guest interaction factor affecting the stability of its clath- 
rates will probably be the diameter of the cavities formed by the host 
molecules. This is consistent with the fact that measured AH p /2AH 
ratios for its clathrate with hydroquinone (1.8) and with water (0.8) are 
approximately inversely proportional to the cavity diameters of the two, 
4.2 and 5.2 A respectively. However, the interactions between methane 
and its host might be enhanced somewhat if the host cavities were lined 
with groups (such as hydroxyl or amino groups) having a high hydrogen 
bonding power as in the case of the apparently enhanced solubility of 
methane in alcohol solvents at a given solubility parameter. In the case 
of clathrates with hydrogen bonding power, one would not have to be con- 
cerned with the solubility parameter of the host. 

1 .0 Natural Gas Hydrates^* 25,26 


Perhaps the most familiar example of clathration to the gas industry are 
the hydrates of methane and other natural gas hydrocarbons. These hydrates 
became of interest to the U. S. gas industry around 1934 when it was noted 
that plugging of natural gas transmission lines was not due to freezing of 
water at 0°F, but to formation of hydrates of hydrocarbon constituents of 
natural gas at temperatures as high as 150°F.^ It was subsequently found 
that gas hydrates exist in two distinct but different crystal structures 
depending on the size of the guest molecule: 

• Structure I hydrates consist of small hydrating (guest) molecules, for 
example, Ar, Kr, Xe, CH^, C^H^, H^S, and CH^Cl. The cubic unit cells 
have 46 water (host) molecules and 8 cavities (voids, cages) where guest 
(hydrating) molecules may be located. The ideal formula (all cavities 
filled) would be — 
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(Solute) *46/8 H 2 0 or (Solute) *5 .75 H 2 0 

• Structure II hydrates consist of larger hydrating (guest) molecules; for 
example, CgHg, CHClg, and C 2 HcC 1. The cubic unit cells have 136 water 
(host) molecules and 8 voids (cavities, cages) for guest molecules. The 
ideal formula (all cavities filled) would be — 

(Solute) *136/8 H 2 0 or (Solute) *17 H 2 0 

The above stoichiometry indicates that methane could be stored as the 
hydrate to the extent of 0.155 g of methane/g water. 

Unfortunately, methane hydrate is not sufficiently stable for purposes of 

storage at ambient conditions. This is shown by Figure D-l which indicates 

that methane hydrate is most stable above 4000 psig at 70°F (27.7 MPa at 21 °C), 

or above 20,000 psig at 100°F^ (138 MPa at 38°C). The equilibrium pressure 

at a given temperature is reduced considerably for the mixed hydrates of 

28 29 

methane with either ethane (Figure D-2) or propane (Figure D-3). Thus, 
the overall vapor pressure of the hydrate at a given temperature can be 
reduced by blending methane with higher hydrocarbons. However, the effect 
with ethane and propane is not sufficient to make the hydrate useful for 
storage under truly ambient conditions. Furthermore, the hydrates of hydro- 
carbons above C 5 are not known. 

Dissolved salts such as sodium chloride significantly decrease the tem- 
perature at which the hydrate is stable at a given pressure (Figure D-4), 
presumably because they interfere with the crystallization (freezing point 
lowering) . 

Thus, it would appear that the mechanical strength of the crystalline 
water host in these hydrates is not sufficient to contain the vapor pressure 
of methane at ambient temperatures. Furthermore, it does not appear that the 
vapor pressure to be contained can be lowered sufficiently by blending the 
methane with higher hydrocarbons or other lower vapor pressure “molecules. One 
wonders then if the strength of the host ice crystal structure can be increased 
in some other way. It was suggested that hydrate forming salts such as MgCl 2 
might have this effect. However, the data on the effect of sodium chloride 
above suggests that the presence of MgCl 2 would also be deleterious. In par- 
ticular, inorganic cations tend to complex with water (i.e., form chemical 
hydrates) through the oxygen atom, thus interfering with the hydrogen bonding 
between hydrogen and oxygen and, thus, with ice formation. On the other hand, 
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Figure D-2. HYDRATE EXISTENCE CONDITIONS FOR METHANE, ETHANE, AND 
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Figure D-3. HYDRATE-FORMING CONDITIONS FOR METHANE-PROPANE 
MIXTURES (From Reference 29) 
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Figure D-4. METHANE HYDRATE PRESSURE-TEMPERATURE CONDITIONS FOR 
VARIOUS SODIUM CHLORIDE BRINE CONCENTRATIONS (From Reference 10) 
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it might be possible to find a suitable inorganic salt which* in the process 
of crystallizing as a hydrate* will clathrate small molecules such as methane. 
However, such a system would not be an extrapolation from the gas hydrate or 
any other known system but rather a new and unexplored area. 

A traditional method of increasing the strength of materials, of course, 
is by reinforcement of the material with fiber (composite materials). A com- 
mon example is glass fiber reinforced epoxy resins. Perhaps then, methane 
hydrate formed in the presence of a matrix of colloidial silica which gels 
into a highly cross-linked fibrous mass could show enhanced vapor pressure- 
temperature relationships. However, this is obviously a speculative sugges- 
tion at the moment. 

Urea and Thiourea Inclusion Compounds* »^ , ^»^ ) ** 

Considerable literature exists on the preparation and properties of urea 
and thiourea clathrates, primarily because of their use in the separation of 
isomers or other mixtures on the basis of difference in molecular size and 
shape. In' particular, straight chain paraffin hydrocarbons (and fatty acids) 
have been separated commercially on the basis of the fact that the urea clath- 
rate structure will accommodate rj-paraffin chains but not branched chaii^ 
hydrocarbons. Thus, a number of excellent reviews are available. 

The urea and thiourea adducts are examples of the tunnel clathrates 
formed by crystallization through hydrogen bonding of the host molecules in 
the form of an extended helix or tunnel around the guest molecule. The tunnel 
diameter is of the order of 0.52 nm for urea and about 0.61 nm for thiourea. 

The adducts are usually prepared basically by crystallization from a suit- 
able solvent (e.g., water or methanol) in the presence of the guest molecule, 
rphexane appears to be the shortest straight chain hydrocarbon that will form 
a urea adduct at ambient conditions. However, urea adducts of propane and 

butane have been prepared by Schlief by lowering the temperature to slightly 

3 1 

below the boiling point of the hydrocarbon with methanol as solvent. It is 
possible that the urea adduct of methane may actually exist and be stable 
under ambient temperatures at elevated pressures. 

Figure D-5 from Schlief^ suggests that if it does exist, the methane/ 
urea adduct would have a composition approximating 1 CH^/3 ureas corresponding 
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Figure D-5. DEPENDENCE OF THE GUEST/HOST RATIO UPON CHAIN LENGTH 
FOR HYDROCARBON/UREA CLATHRATES (From Reference 31) 



to a weight ratio of 0.089 g/g which is of the same order of magnitude as that 
for adsorption on activated carbon at 3.6 MPa. 

2.0 Phenolic Clathrates 1 ,2 » £ > 12 9 ^ 

Of considerable interest is the clathration behavior of phenols in general, 
and hydroquinone in particular. In the presence of small gaseous molecules 
such as methane (or O 2 , N 2 , C 2 H 2 , ^ 3 ^, HC1, SO 2 , Ar, Kr, Xe), hydroquinone 
crystallizes in the metastable "g" modification in a ratio of 1 cavity to 
every 3 hydroquinone molecules. This crystal modification, however, is stable 
only if a certain fraction of the cavities are filled with guest molecules 
which are neither too small (He) nor too large (CCI 4 ). The cavity has a diam- 
eter of 3. 95 A when occupied by small molecules comparable to 3. 8 A for methane. 
However, larger molecules can be accommodated by distortion of the cavities. 
Although the crystal structure is held together by hydrogen bonding in a man- 
ner similar to ice, the stability of the hydroquinone clathrate is far great- 
er: the CH^ * 3 C^H^(OH )2 complex is quite stable at ambient temperature and 

pressure. X-ray diffraction studies indicate that the clathrate involves the 
formation of the typical hexagotially hydrogen bonded structure formed from six 
hydroxyl groups as discussed above. 

The clathrate can be prepared by crystallization from a saturated solu- 

32 

tion of hydroquinone in ethanol in the presence of methane pressure or from 

34 

the vapor phase comprising hydroquinone vaporized in a methane carrier gas. _ 

In both cases, the fraction of the cavities filled with methane is highly 
dependent on the methane pressure (Figure D- 6 ); thus, complete filling to 
yield the one to three stoichiometry requires a pressure of about 100 atm 
(10.4 MPa). The methane can be recovered for use by dissolution in water or 
ethanol or by heating (melting). No data was found on the thermal stability 
of the clathrate. However, the melting point of the modification of hydro- 
quinone is 50°C (123°F). Thus, it is probable that it will decompose at about 
that temperature. 

Unfortunately, the gravimetric capacity of the methane/hydroquinone clath- 
rate is only 0.0485 g CH^/g hydroquinone due to the high molecular weight of 
the host. 

Phenol and several substituted phenols also form clathrates with methane 

1 3 A 

and other gases, including rare gases, ethane, propane, and the butanes. * 
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PRESSURE, atm 


Figure D-6. FRACTION OF CAVITIES FILLED WITH METHANE 
AS A FUNCTION OF PRESSURE (from Reference 32) 
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These clathrates are also presumed to involve the typical hexagonal hydrogen 
bonded oxygen structure from six phenol groups together with the typical 1/3 
stoichiometry. However, the data for methane clathrates with the phenols pre- 
pared at -196°C indicate a stoichiometry closer to 1 CH^/2.4-2.6 phenol 
molecules (instead of 1/3 as with hydroquinone) • This may indicate that more 
than 1 methane molecule can be accommodated in the cage structure of the less 
lightly bound phenol group. This hypothesis is supported by the fact that 
hydrocarbons up to n-butane can also be clathrated by the phenols although at 
levels less than 1 molecule per cavity. ^ 

The data of Lahr and Williams on rare gas /phenol clathrates indicates 

that the CH^/phenol clathrate would be much less stable than the hydroquinone 

35 

analog at ambient temperature. 

Dianin's Compound^ > 12-14 

Dianin’s compound is the product of the condensation of mesityl 

oxide with phenol and is another example of the hexa-host class of clathrates 
described above, based on the hydrogen bonding of the phenolic hydroxyl groups. 
The cavities formed are hour-glass shaped, about 11A in length, 4. 2 A at the 
waist, and 6. 4 A at the two widest points. Although best known for inclusion 
of solvent molecules, it also forms clathrates with fixed gases including 
methane (Barrer and Shanson).^ Dianin’s compound differs from other phenolic 
clathrates in that its g clathrate structure is stable whether the cavities 
are filled or not. Normally, three mole of phenolic or other hexa-host type 
are associated with each cavity. The apparent saturation capacity of Dianin’s 
compound for methane is 6 CH^/6 Dianin compounds or 0.06 g/g which may indi- 
cate that three CH^ molecules occupy each hour shaped cavity. However, this 
was measured at the boiling point, 182. 5 °C, and Barrer indicates that the 
apparent saturation capacity varies with temperature (at atmospheric pres- 
sure). ^ This suggests that such capacities might also exist at ambient 
temperature with elevated pressure. 

Even more important, the data indicate that several CH^ molecules can 
exist in a clathrate cavity, at least in more or less cylindrical cavities 
having diameters not too much larger than the kinetic diameter of methane. 
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Cy clodext r ins * > 2 , 9-1 3 

The cyclodextrins are cyclic oligosaccharides containing six (a), seven 
(ft), or eight (v) glucose units. The molecules in solution are presumably 
doughnut shaped with hole diameters of 6 A, 8 A, or 10-11 A, respectively, for 
the a, g, and v forms. The voids in the center have been described as having 
a high electron density (Shaeffer and Dorsey),*® capable of complexing with 
various molecular species since no cage structure is evident. Furthermore, 
clathration has been shown to change the oxidation-reduction potential of 
certain molecules such as methylene blue.*® 

Alpha cyclodextrin (6 A) also forms clathrates with various fixed gases, 
including methane at ambient temperature-pressure to the extent of 0.014 g 
CH^/g-ot cyclodextrin. However, these clathrates are stable only in the dry 
state and, therefore, may involve cage structures. Nonetheless, this illus- 
trates that cages as large as 6A (about 2A greater than the kinetic diameter 
of methane) can exhibit relatively high stability at ambient conditions, prob- 
ably due to the inherent stability of the "beta” structure. The latter is 
identical to the "alpha” structure, at least in the dry state, since the cav- 
ity is held together by covalent bonds rather than hydrogen bonds. However, 
the vaporization of methane from the clathrate may also exhibit an enthalpy 
change greater than AH vap in spite of the minimal coordination number due to 
the high electron density of the oxygen atoms lining the cavity walls. 

The six, seven, and eight raerabered cyclodextrins are natural products 
produced by enzymatic action on amylose. However, one wonders whether a five 
membered analog could be produced synthetically. If so, the hole diameter 
would be of the order of 4A (by extrapolation), just right for inclusion of 
methane (d = 3. 8-4.1 A). Furthermore, it may be possible to synthesize cyclic 
polyethers (from ethylene oxide) tailored for inclusion of methane. 

A further possible approach would involve the complete synthesis or modi- 
fication of existing cyclic polyethers or cyclodextrins with low solubility 
parameter detergent-like tails to produce inclusion hosts soluble in solvents 
which also dissolve methane. In this manner, the solvent capacity for methane 
might be enhanced by the action of the dissolved inclusion hosts. 
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3 .0 Encapsulation in Preformed Molecular Sieves^ > ^ ^ 

In the broadest sense, zeolites, molecular sieves, and activated carbons 
can be considered analogous to the host crystal of a clathrate except that the 
structure is preformed and stable without the included guest molecules. Such 
molecules have high internal surface areas contained within cages, or pores 
(tunnels) having precisely sized entries. Such materials are used not only 
for adsorption but also for separation of molecules having different kinetic 
diameters. Under ordinary conditions, the guest molecules must be smaller 
than the kinetic diameter of the host; otherwise, they will be excluded. It 
has recently been found that such molecular sieve structures can also be used 
for encapsulation or clathration of guest molecules larger than the charac- 

07 

teristic pore opening size. This is accomplished at high pressures and 
elevated temperatures such that the molecules are forced into the slightly 
expanded pore openings. Once inside, and after the temperature has been 
quenched to ambient levels, the molecules are trapped at high pressure inside 
the sieve. Such a system tends to be quite stable at ambient temperature and 
pressure until reheated, destroyed by acid, or exposed to a stronger absorbate 

9 U2O) • 

More specifically, it has been found that inert gases such as argon, 
krypton, and methane can be stored effectively in K-A-type zeolites at 25°C 
after encapsulation at 300° to 400°C and 200 to 400 MPa (2050 to 4100 atm). 
Data for encapsulation of CH^ and ethylene as well as krypton and argon are 
shown in the following table. 

Table 1. ENCAPSULATION OF GASES IN ZEOLITE A 17 


Gas 

a (A) 

Temperature, 

(°C) 

Pressure, 

(MPa) 

3 

cm 

STP/s 

After t 
Days at 
25°C 

ch 4 

3.8 

350 

268 

105, 

t = 4 

98, t = 37 

c 2 h 4 

3.9 

250 

84 

81, 

t = 2 

76, t = 35 

Ar 

3.4 

350 

268 

109, 

t = 3 

77, t = 79 

Kr 

3.6 

350 

436 

90. 

5, t = 1 

90.5, t = 30 
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Thus, at 350°C and 270 MPa, 105 cc STP/g are encapsulated. This corresponds 
to 0.075 g/g, which is of the same order of magnitude as for adsorption of GH^ 
on zeolites at 3.6 MPa and 25°C. It should also be noted that after 37 days 
at 25 °C and 1 atm, less than 7% of the CH^ is lost. However, because of the 
drastic conditions required for the encapsulation of CH^ into zeolite A, the 
process is obviously not feasible. Furthermore, the temperature required to 
decapsulate is undoubtedly too high (350°C). 

Figure D-7 shows the pore sizes of a number of zeolites relative to the 
kinetic diameters of a number of gases including CH^. It will be noted that 
zeolite KA has a pore size opening of about 3A compared to 3 *8 A for methane, 
while LiA at 3. 3 A or NaA at 3 ,5 A are closer to that of methane. In addition, 
other molecular sieve materials, including activated carbons and the newly 
developed alumino phosphates, could be of interest The average pore 
size of an activated charcoal can be varied over wide limits, depending not 
only on the conditions of activation, but particularly on the basic structure 
of the char substrate before activation. For example, Saran charcoals made by 
pyrolysis of polyvinylidene chloride have been shown to have a uniform slot- 
like pore structure between graphite-like layers of carbon. The slot-like 
pores are about 10 A in length and thin enough to restrict but not exclude the 
flow of branched chain hydrocarbons (neopentane) as compared to planar (ben- 
zene) or linear (rrpentane) hydrocarbons. A molecular sieve carbon made by 
Takeda Chemical Industries, Osaka, Japan, exhibits an even smaller uniform 
pore structure of only 5A.^ It is interesting to note that the heats of 
adsorption of various gases on these materials were twice those measured for 

/ Q 

graphitized carbon black. The increased AH (adsorption) was attributed to 
the fact that adsorbed species interact with carbon layer planes on both 
sides, while molecules adsorbed on surfaces interact with only one carbon 
layer plane. Thus, even though the pore openings of these carbons are too 
large for encapsulates, the small pores offer the advantage of higher heats of 
adsorption which will result in the achievement of a given surface coverage at 
a lower pressure. 

A new class of molecular sieves, the alumino phosphates, have recently 
been introduced by Union Carbide Corp.^*^* None of t^iese materials, which 
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have so far been described, is particularly suited for encapsulating methane. 
Nonetheless, it seems a reasonable possibility that alumino phosphate or even 
zeolite molecular sieves could be tailored specifically for that possibility. 


Cross-Linked Sorbents Formed in the Presence of Templates 


The addition of templates or inclusion bodies during formation to control 

the size of the cavities and, therefore, the specificity of various adsorbents 
has been done with a number of materials including: 




Zeolites 

involved 1 


nd alumino phosphates in which the templates are the cations 

, 


• Silica gels in which the gels are gelled around the template molecules, 
which are usually relatively larger, water-soluble molecules such as 
butyl orange dye 1 ^ 

• Modified dextrans which are macromolecular polysaccharides cross-linking 
to varying degrees in the presence of a suitable solvent and are used for 
gel permeation separation of molecules on the basis of molecular size. 

Thus, it may be possible to develop adsorbents or absorbents tailored 
specifically for methane by controlled cross-linking of organic as well as 
inorganic polymers in the presence of high pressures of methane (or other 
templates of a similar size, including propane and butane). Examples would 
include: 

• Polymerization of isobutylene with suitable cross-linking agent (e.g., 
diisobutylene) in a propane or butane solvent 

• Three-dimensional polymerization of a silicon monomer or oligomer with a 
suitable cross-linking agent 

• Cross-linking of various polymers including a hydrogen bonding polymer 
such as cellulose (with its appended hydroxyl groups) or polyacrylamide 
(with its amido or substituted amido groups) or polymers having a solvent 
affinity for methane such as polyisobutylene swollen by a low molecular 
weight solvent which in the dry state should contain a multiplicity of 
cavities lined with either hydrogen bonding groups or methyl groups. 

• Vapor phase formation of cross-linked silica gels in the presence of high 
concentrations or pressures of methane. This might be accomplished by 
saturating a methane stream with a suitable silicon compound (e.g., SiCl^ 
or Si(O-CH^)^) followed by vapor phase hydrolysis or vapor phase thermal 
decomposition. 
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6.0 Other Clathrate Systems 1 > 2 ^ ** 

There are many other molecules which have been shown to clathrate guest 
molecules including: 

• Phosphonitrites 

9 Werner complexes 

• Other metal ion-ligand combinations 

• Steroids 

• Cholesterol 

9 Hexamethyl isocyanide chloride 

• Methyl naphthalene 

• Cycloveratril. 

These clathrate types will not be described here since no methane clath- 

rates were attempted, even though higher molecular weight guests were included. 
However, such clathrate hosts should not be excluded from any further studies 
with methane since the latter was never even tried. For example, it is pos- 
sible that some of these would crystallize differently in the presence of 
methane or may contain layer-like cavities which would accommodate methane as 
well as larger guest molecules. In any case, it may be possible to vary the 
molecular structures of the host to give the methane clathrate. Of particular 
interest are the metal-ion-ligand combinations since they are capable of wide 
variation in both the ion and the ligand used. 

7 .0 Summary 

/ 

We have seen that methane clathrates exist which separately exhibit the 
various properties required for on-board storage, although not together in the 
same clathrate. Furthermore, a study of the literature suggests that nearly 
any chemical species near or below its freezing point has the potential of 
acting as a clathrate host under the right conditions and in the presence of a 
reasonable concentration of a suitable guest molecule. Thus, it is possible 
that a methane clathrate having all of the properties required for on-board 
storage can be found or developed although this may entail a compromise be- 
tween stability, capacity, and methane release properties. 
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We also recognize that clathrates have properties akin to solutions on 
one hand and to molecular sieve-absorbate interactions on the other. This is 
illustrated by a comparison of intracrystalline free volumes available for in- 
clusion of g-clathrate hosts with the volumes available in zeolites as follows 

HydroquinOne /argon ~0.05 cc/cc 

Urea/paraffins 0.37 

Thiourea/hydrocarbons 0.41 

Gas hydrates 0.46 

Zeolites 0 .18-0.54 

Thus, there is a large overlap in the potential inclusion volumes of zeolites 
and clathrates which is comparable to that in solutions. Inclusion of up to 
about 50% by volume are potentially possible by all three approaches, although 
the rules governing the accomplishment thereof are different. 

For solutions, the primary influence appears to be the solubility param- 
eter which appears to be inversely proportional to molecular weight so that 
50% (by volume) solutions can only be achieved using solvents with molecular 
weights close to that of methane itself. This is not feasible for our pur- 
poses. With zeolites (and other adsorbents) for which the host structure is 
inherently stable by itself, the inclusion is limited at a reasonable pres- 
sure, e.g., 3.6 MPa, by the heat of adsorption of methane on the surface to 
less than a monomolecular coverage of the surface, which is usually less than 
potentially possible with inclusion compounds. The latter, however, can only 
be accomplished at high pressures or with large molecules at lower pressures. 

With clathrates, stability of the clathrate depends on the stability of 
the 3 -crystal lattice relative to that of the usual orform and upon the magni- 
tude of the heat of interaction between host and guest. Both contribute to 
the difficulty of removal of the guests. With most of the methane clathrates 
so far found, the potential saturation of the free volume available approaches 
100 %, probably because the clathrate is not discovered unless stability is 
achieved. As indicated above, stability under reasonable conditions is usu- 
ally not achieved unless the cavity dimensions approach that of the guest 
molecule so that a relatively high coordination number for interaction is 
achieved. 
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Suggested Techniques for Future Clathration Host Development: 

Although there is considerable work in progress toward developing the 
know-how for tailoring new clathrate hosts for specific guest molecules, most 
of the work is oriented toward biological and catalytic applications using 
enzymes and ionic species. No references have been found indicating ongoing 
work directed toward tailoring clathrate host molecules for small and rela- 
tively inert molecules such as methane or the inert gases. However, should 
such a program be undertaken, the following precepts are available to guide 
such development: 

• Clathrate formation depends on the ability of the host molecules to form 
a crystal habit which has cavities large enough to accommodate the speci- 
fied guest molecules. This is usually, although not necessarily, a crystal 
habit other than the normal crystal structure formed in the absence of 
the host. This ability would be difficult to predict, a priori, although 
it might be possible for a well versed crystallographer to develop such a 
predictability on the basis of molecular structure. 

Indeed, this is essentially the approach used by MacNicol, et al. ,^ who 
noted that many of the known clathrate structures involve the formation of a 
hexagonal ring of hydrogen-bonded oxygen atoms from six phenolic hydroxyl 
groups as shown below. 





Structure A 

2 

The clathrate hosts which show this behavior include Dianin’s compound, re- 
lated compounds, hydroquinone , ^ phenols,^ and possibly water. MacNicol, et 
al., then demonstrated that the clathration ability of such hexagonal struc- 
tures could be simulated using the permanent hexagonal structure of certain 
hexa substituted benzenes: 
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Structure B 

where R-Z-radicals included a variety of substituted oxygen and sulfur radi- 
cals. However, the guest molecules bound by these clathrate hosts are larger 
than methane, indicating relatively large cavities. 

By the same token, new clathrate hosts might also be derived by varying 
the nature of R- in Structure A to include substituted phenol and also various 
aliphatic, alicyclic, and olefinic species, substituted or unsubstituted. 
Furthermore, the formation of suitable hydrogen-bonded hexagonal structures 
can also be formed with elements other than oxygen from Rows 5A and 6A of the 
periodic table, in particular through nitrogen atoms in amides. In addition, 
the formation of a hexagonally arrayed hydrogen-bonded structure is not an 
exclusive prerequisite, since urea forms tunnel clathrates with straight chain 
paraffins by forming a hydrogen-bonded two-dimensional helix spiral. (tunnel) 
around the paraffin molecule. Thus, the primary prerequisite is the hydrogen- 
bonding (or other complex forming) ability of the host while the secondary 
prerequisite is a molecular geometry (at least partially unpredictable) which 
will crystallize in the presence of a guest to a structure with cavities sized 
to accommodate the guest. 

Methods for Increasing the Stability of Clathrates and Improving Methane 
Release Properties: 

Stability is of major importance since it controls not only whether the 
methane can be contained at reasonable pressures under ambient conditions but 


D-24 



also the conditions under which the methane can be released. The hydrate of 
methane for example is not stable enough since it requires a pressure of 28 
MPa to exist under ambient temperature conditions. On the other hand, the 
hydroquinone /me thane clathrate is too stable. Thus, although it can be stored 
at ambient temperature and pressure, it requires an elevated temperature 
(~50°C) or the introduction of a solvent to release the methane. Obviously, a 
compromise is in order so that the clathrate will be stable at ambient condi- 
tions at a reasonable pressure of 1.4 to 3.6 MPa, in order that methane re- 
lease can be achieved and controlled by pressure reduction. 

The parameters controlling stability would include: 

• Hydrogen bonding power of the groups through which the crystallization 
occurs, as determined by the nature of the groups themselves, as well as 
the size and electrophilic character of the other groups in the molecule. 
Of particular interest in this regard would be the amides and substituted 
amides. 

• Geometry and symmetry of the host molecule as it affects the structure of 
the crystal habit formed. The importance of symmetry is suggested by the 
fact that the hydroquinone (p-dihydroxybenzene) gives a highly stable 
clathrate with HBr, whereas the meta analogue does not exist at atmo- 
spheric pressure. 

Suggestions for Increasing Storage Capacity of Methane: 

The capacity of a clathrate for storing methane depends on two factors: 

• How many host molecules are required to provide one guest "cage” 

(i.e., the unit cell) 

• Molecular weight of the host molecule. 

The structure of the unit cell and, therefore, the number of host molecules 
per guest molecule is dependent in a very complex way on various crystallo- 
graphic factors and is beyond prediction at this time, except perhaps for the 
"hexamer" clathrates (described above) which tend to give a 3/1 host to guest 
ratio. However, the unit cell is likely to contain a reasonably small number 
of host molecules, e.g., 

• 3 in the hydroquinone-methane clathrate and other ”hexamer ,, hosts 

• 5.75 in the methane hydrate. 

Thus, molecular weight may be the factor of greatest importance in determining 
capacity. 
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On the other hand, it may be possible that some higher molecular weight 
molecules with complex symmetries may form clathrate structures in which the 
host/guest ratio is significantly less than the 3/1 found for the hexamer 
types, or that form structures with cavities sufficiently large to accommodate 
two or more methane molecules, as in the case of the urea tunnel adducts. In 
the latter case, diffusion would be the primary barrier to decomposition and 
the system would tunction in a manner analogous to that of absorption in a 
zeolite. Most of the clathration entities forming larger cavities (and there 
are a lot of them, including Dianin’s compound, the Werner complexes, deoxy- 
cholic acid, and urea and thiourea, etc.) have been evaluated only with larger 
guest molecules (usually the solvent used) at atmospheric pressure and not 
with permanent gases such as methane under elevated pressures. For example, 
urea adduction of the homologous n-paraffin series has been extended downward 
to include propane and butane, which are stable at atmospheric pressure only 
at subzero temperatures. However, it is possible that adduction could be 
extended to include methane and ethane as well at ambient temperatures and 
elevated pressures. 

8. Future Research Recommendations 


It is apparent from the above discussion that any future work toward the 

development of clathrate systems for storage of natural gas will be highly 

empirical in approach. Such a program would involve: 

1) Evaluation of the formation and stability of clathrates of the known host 
systems toward methane at pressures up to 7.0 MPa. In particular, this 
would include the urea and thiourea systems and those whose cavities are 
presumed to be too large for methane as suggested above (particularly the 
ligand cross-linked metal ion, e.g., Werner complexes). 

2) Evaluation of the effect of variations in structure on stability and 
capacity. Of particular interest would be: 

• N-substituted ureas and thioureas to enhance hydrogen bonding strength 
and, therefore, stability of the adducts 

• Diphenols including the effect of ring substitution on stability of 
hydroquinone clathrates and the use of meta and ortho diphenols 

e Variations based on cyclodextrin, including synthesis and evaluation of 
the five membered oligomer and conversion of the inwardly-oriented 
hydroxyl groups to ethers of varying size, amides, etc. to decrease the 
hole size, and/or the interaction chemistry with methane. 
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3) Determine the effect of more drastic changes in the molecular structure 
of amides on their clathration behavior, e.g., 

0 0 0 0 0 0 

I I I I I I 

NH 2 -C-C-NH, NH 2 -C-(CH 2 ) n -C-NH 2 , CH 3 -C-NH 2> CH 3 (CH 2 ) 2 -C-NH 2 

0 0 0 0 0 

I II II 

C, H-C-NH , NH -C-NH-C-NH 0 , NH 0 -C-C-NH„ 

0 0 l L l L l 


4) Synthesis and evaluation of cyclic polyamides (analogous to the crown 
ethers but with greater hydrogen bonding propensity) 

5) Evaluation of systems based on aliphatic alcohol groups such as glycols 
and substituted glycols, glycols and substituted ethers and polyethers 

6) Evaluation of systems based on aliphatic acids such as oxalic acid, 
acetic acid, chloroacetic acid, malonic acid, chlorinated malonic acid, 
maleic acid, and succinic acid derivatives. 


D-27 



REFERENCES FOR APPENDIX D 


1. Bhatnagar, V. M., Clathrate Compounds, 1970, Chemical Publishing Co., NY 

2. Makin, E. C., "Clathration, " Encyclopedia of Chemical Technology, Vol. 6, 
p. 179-189, Kirk and Othmer, Eds., New York: Wiley (1978). 

3. Davidson, D. W. , Water: A Comprehensive Treaties, Vol. 2, Plenum Press, 
New York (1973). 

4. Hagan, Sister Martinette, Clathrate Inclusion Compounds, 1962, Reinhold 
Publishing Co., NY 

5. Powell, H. M. , "Clathrates", Chapter 7 in Non-Stoichiometric Compounds . 
Mandelcorn, L., Ed., Academic Press, NY (1964). 

6. Fetterly, L. C., "Organic Adducts," Chapter 8, Ibid. 

7. Barrer, R. M., "Inorganic Inclusion Complexes," Chapter 6, Ibid. 

8. Senti, F. R. and Erlander, S. R. , "Carbohydrates," Chapter 9, Ibid. 

9. Staveley, L. A. K., "Physics and Chemistry of Inclusion Complexes," Chap- 
ter 10, Ibid. 

10. Schaeffer, W. D. and Dorsey, W. S., "Clathrates and Clathrate Separa- 
tions,” Chapter 3, Advances in Petroleum Char and Refining VI, 119-167 
(1962). 

11. Cramer, F. D., "Inclusion Compounds," Rev, of Pure and Applied Chem. , 

Vol. 5, 143-164, (1955). 

12. Mandelcorn, L., "Clathrates." Chem. Rev. 59, 827-839 (1959). 

13. Frank, S. G., "Inclusion Compounds," J. Pharmaceutical Sciences, 64, #10, 
1585-1604, (1975). 

14. MacNicol, D. D. , McKendrick, J. J., and Wilson, D. R., "Clathrates and 
Molecular Inclusion Phenomena," Chemical Soc. Reviews, 65-87 (1978). 

15. Child, W. C., "Molecular Interactions in Clathrates: A Comparison with 
Other Condensed Phases," Quart. Rev., Vol. 18, p 321-346. 

\ 

16. Van der Waals, J. H. and Platteeuw, J. C., "Clathrate Solutions" in 
Advances in Chemical Physics. Vol. II, Edited by Prigogine, I. Inter- 
science Publishers, Inc., NY, 1959, p. 1-57. 

17. Breck, D. W., "Crystalline Molecular Sieves," J. Chem. Ed. 41, 678-689, 
(1964). 


D-28 



18. Breck, D. W. , Zeolite Molecular Sieves, John Wiley & Sons: New York, NY, 
1974. 


19. Cram, D. J., "Cavitands: Organic Hosts with Enforced Cavities," Science 
219, 1177-1183, (1983). 

20. Cram, D. J. and Cram, J. M. , "Host-Guest Chemistry," Science, 183, 803- 
809, (1974). 

21. MacNicol, D. D. , and Wilson, D. R. , "New Strategy for the Design of In- 
clusion Compounds: Discovery of the "Hexa-hosts" ," J. Chem. Soc. Comm. p. 
494-495, (1976). 

22. Hardy, A. D. U., MacNicol, D. D., and Wilson, D. R», "A New Approach for 
the Design of Inclusion Compounds," J. Chem. Soc. Perkin II, 1011-1019, 
(1970). 

23. MacNicol, D. D. , Hardy, A. D. U., and Wilson, D. R. , "Crystal and Molecu- 
lar Structure of a "Hexa-host" Inclusion Compound," Nature, 266, 611, 
(1976). 

24. Scott, M. I., Randolph, P. L. , and Pangborn, J. B., "Assessment of 
Methane Hydrates." Final report for period Dec. 1978 through June 1980, 
GRI-79/0070, Contract #5011-310-0097, Oct. 1980. 

25. Hammer schmidt, E. G., Knapp, K. R. and Perkins, C. L., "Gas Hydrates and 
Gas Dehydration," in Gas Engineers Handbook . New York: Industrial Press, 
1969. 

26. Parent, J. D. , "The Storage of Natural Gas as Hydrate," IGT Res. Bull. 

No. If Chicago, 1948. 

27. Hand, J. H. , Katz, D. L., and Verma, V. K. , "Review of Gas Hydrates With 
Implication for Ocean Sediments," in I. R. Kaplan, Ed., Natural Gas in 
Marine Sediments, 179-94. New York: Plenum Press, 1974. 

28. Baker, P. E., "Experiments on Hydrocarbon Gas Hydrates in Unconsolidated 
Sand," in I. R. Kaplan, Ed., Natural Gas in Marine Sediments , 227-34. 

New York: Plenum Press, 1974. 

29. Katz, D. L. et . al. , "Water-Hydrocarbon Systems," in Handbook of Natural 

Gas Engineering. New York: McGraw-Hill, 1959. 

30. Kobayashi, R. et . al . , "Gas Hydrate Formation with Brine and Ethanol 
Solutions," in Proceedings Thirtieth Annual Convention, Natural Gasoline 
Association of America, 27-31, Tulsa, Okla. , April 25-27, 1951. 

31. Schlief, V. H., "Harnstaf fein Schlusverb in dungen von q-Butan und 
v-Propan," J. Prakt. Chem. , _4_, 335-339 (1955). 

32. Peyronel, G. and Barbieri, G., "On Some New Clathrates of Hydroquinon, ” 

J. Inorg. and Nuclear Chem. , 8, 582-585 (1958). 


D-29 



33 


. Coutant, R. W. , "Solventless Preparation of Hydroquinon Clathrates , " J. 
Org. Chem. , 39, 1593-1594 (1974). 

34. Barrer, R. M. and Shanson, V. H. , "Clathration by Para-Substituted 
Phenols,” J. Chem. Soc., Faraday, T. t 2348-2354, (1976). 

35. Lakr, P. H. and Williams, H. L., "Properties of Some Rare Gas Clathrate 
Compounds," J. Phys. Chem., 63, 1432-1433, (1959). 

36. Barrer, R. M. and Shanson, V. H. , "Dianin's Compound as a Zeolitic 
Sorbent," j. Chem. Soc. Comm., 333-334, (1976). 

37. Sesny, W. J. and Shaffer, L. H., U.S. Patent No. 3,316,691 (1967). 

38. Fraenkel, D., "Encapsulate Hydrogen," Chemtech., 11, 60-65 (Jan. 1981). 

39. Dacey, J. R. and Thomas, D. G. , "Adsorption on Saran Charcoal," Trans. 
Faraday Soc., 50, 740-748 (1954). 

40. Wilson, S. T., et . al. , "Alumino Phosphate Molecular Sieves: A New Class 

of Microporous Crystalline Inorganic Solids," JACS, 104, 1146-47, (1982). 

41. Haggin, J., "Alumino Phosphates Broaden Shape Selective Catalyst Types," 
C&EN, June 20, 1983, 36-37. 

42. Kawazoe, K. , et . al. , "Correlation of Adsorption Equilibrium Data of 
Various Gases and Vapors on Molecular-Sieving Carbon," J. Chem. Eng, of 
Japan, Vol. 7, p. 158-162 (Eng.), (1974). 

43. Kazuyuki, C., Suzuki, M., and Kawazoe, K. , "Adsorption Rate on Molecular 
Sieving Carbon by Chromatography," A.I.Ch.E. Journal, Vol. 24, March, 
237-246, (1978). 


D-30 



APPENDIX E 


EVALUATION OF DISSOLUTION AS A MEANS OF 
STORING NATURAL GAS 

SOLUTIONS 


1 • Background 

The concept of an ideal solution is used to describe the behavior of 
actual solutions as a first approximation, similar to the manner in which the 
concept of the ideal gas is used to describe real gases. In both cases, the 
degree of accuracy of the approximations are reasonable for many chemical 
entities, at least within certain limits, i.e., for dilute solutions or for 
gases at low pressures.* "In both cases, the idealized limiting behavior can 
be either: 1) defined thermodynamically by means of empirical expressions, or 
2) derived from idealized models of molecular systems."* The thermodynamic 
definition of an ideal solution is: "one in which the activity equals the mole 
fraction over the entire composition range and over a non-zero range of temper- 
ature and pressure." Thus, for solutions of gases which approximate the ideal 
gas laws reasonably closely, the ideal solution is defined by Raoult’s law: 


- p°X 
1 11 

and 


P° X 
z *2 


( 1 ) 


where: 

i i 

P£ or P£ = Ideal partial vapor pressures of Components 1 and 2 

Xj or X 2 = Mole fraction of Components 1 or 2 in liquid mixture 

?1 or ^2 ~ Equilibrium vapor of Components 1 and 2 over pure materials. 

Raoult’s law, in this case, was deduced on the assumption of zero heat of mixing 
and an entropy of mixing that is independent of temperature and of peculiarities 
of molecular sizes and shapes. Thus, an ideal solution is defined in terras of 
zero heat of mixing. A corollary of this definition is Henry’s law which states 
that the equilibrium value of the mole fraction of a gas dissolved in a liquid 
is directly proportional to the activity (or partial pressure if the gas approx- 
imates ideality) of that gas above the liquid surface. 

Xj - K H P 1 (2) 
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Where: 


Kjj = Henry’s law constant. 

Henry’s law is not restricted to ideal solutions, at least if corrections for 
deviations of the gas from ideality are taken into account. Thus, solubil- 
ities of gases in various solutes are often defined in terms of the Henry’s 
law constant determined at 1 atmosphere pressure (or partial pressure). 

According to Raoult’s law, methane can be stored at ambient temperatures 
at substantially reduced pressures as an ideal solution in a liquid of low 
volatility. However, the vapor pressure of liquid methane at 25°C (extrapo- 
lated) is 29.3 MPa (4248 psia).* Thus, methane stored at 25°C as a 50 mole 
percent solution in a non-volatile solvent will still exert a pressure in 
excess of 14 MPa (2000 psi). To reduce the ambient storage pressure to 3.6 
MPa at 25 °C requires reducing the mole fraction of methane in the solution to 
somewhat less than 0.135; reduction to 1.8 MPa, to less than 0.07 mole frac- 
tion. It is apparent, therefore, that the utility of storage as an ideal 
solution is not a particularly promising approach, unless a noriideal solvent 
can be found such that the solution with methane exhibits a large negative 
deviation from Raoult’s law. 

2. Regular Solution Theory* 

Many non-ideal solutions have sufficient thermal energy to overcome the 
tendency to segregate, thus displaying nearly ideal entropy of mixing. Such 
solutions have been designated by Hildebrand and Scott as "regular" solu- 
tions.* Differences in solubility (i.e., differences in the change of free 
energy of mixing, AF^) compared to ideality will depend on differences in the 
heat of mixing, AH^. On the basis of regular solution theory, Hildebrand 
developed an equation for the solubility of gases in the form: 

V 

log X 2 = log X 2 - 2.303 RT ^ S 1 S 2^ ^ 

3 

Gjaldbaek and Hildebrand subsequently modified the equation for Flory-Huggins 
mixing as follows: 
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log X 2 log X 2 - 2 ^ 303 RT 


V V 

(Sj- S 2 ) 2 - [log ^ + 0.434(1 - ^)[ 


(4) 


Where: 


X 2 = 
Vi,V 2 = 
R 
T 

s l> s 2 “ 


Equilibrium mole fraction of Component 2 in a liquid mixture 

Ideal solubility of Component 2 in mixture 

Molar volumes of components 

Gas constant = 1*9865 

Temperature of solution, °K 

Solubility parameters of Components 1 and 2. 


The solubility parameter, also called the cohesive energy density, S, in turn 
is defined as follows: 


AH. V - RT 1/2 

^ ) 


(5) 


Where: 

AHj = Heat of vaporization of Component 1 at solution temperature 
R,T,Vj = As above. 

The ideal solubility, x|, can be estimated by the following equation: 

l° g X 2 = 2.303 R ^T “ T7^ ^ 

B 

Where: 

T = Solution temperature, °K 

Tg = Boiling point, °K. 

Thus, the primary parameters controlling the solubility of a non-polar gas 
such as methane in non-polar solutions is presumably the temperature of the 
solution relative to the boiling point of the gas, the molal volumes of the 
components and, in particular, the solubility parameter, S, as defined above. 
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It will be noted that if roe thane (or any other non-polar gas) is dissolved in 
a solvent having the same molal volume and the same solubility parameter as 
methane , the last two terms in the equation would be zero and the solubility 
would be that calculated for an ideal solution. 

The Hildebrand treatment of regular solutions has been modified empiric- 
ally to yield better correlations with non-associated polar solvents by Yen 
and McKetta.^ However, the differences in actual solubility covered by these 
correlations are marginal compared to the effects sought in this review. 

Other reviews of the literature on gas solubility in liquids include Battino, 
Lawson, ^ and Pierotti.^ 

Considerable data on the solubility parameters of various molecular 
species are available in the literature, most of which are based on heats of 

vaporization measured at 25°C.* >^,5,6,9 The data presented by Hildebrand and 

12 V 

Scott y specifies the temperature at which aH was measured. However, the 

temperature at which the AH V were derived is not necessarily clear in other 

sources although 25°C is usually implied. In particular, sufficient data for 

methane and ethane (aH V + V) were not found to give good estimates of S at 

25°C. Furthermore, in some literature sources, solubility parameters were 

estimated from solubility data in solvents for which S values are avail- 
12 7 8 

able. * * * However, in this case, variations attributed to S may actually 
be the result of variations in molar volume (i.e., in the Flory Huggins cor- 
rection) • 

Alternate methods of estimating the solubility parameter are also avail- 

12 5 

able as outlined by Hildebrand and Scott * and by others. 

o 

According to the data cited by Hildebrand and Scott, methane behaves in 
most solutions like other non-polar molecules, subject primarily to the van 
der Waals forces. For example, the solubility of methane in various solvents 
correlates nicely with that of other non-polar gases, including H 2 , N 2 , CO, 

°2* C0 2» anc * the rare gases, on the basis of their "force constants,” the 
parameters of the intermolecular energy function as expressed by Lennard- 
Jones. Furthermore, the solubility of methane correlates rather well with the 
solubility parameter of the solvent, which is a measure of the cohesive energy 
density of the solvent. Exceptions *to this, which tend to indicate solvation 
or some degree of polar or acid-base interaction (as in the case of CO 2 " 
benzene mixtures ) 9 have not been found for methane solutions. 
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The efficacy of Equation 4 for estimating the solubilities of a number of 

non-polar gases, including CH^, in various solvents has been tested by Gjaldbaek 

3 7 8 

and associates with reasonably positive results, * * However, these results 
were dependent on estimations of the solubility parameter of methane (and other 
gases near or above their critical points) from Equation 4 using solubility 
data in other solvents of known solubility parameters. Their estimate of S 
for methane on this basis is 6.2, far greater than our own estimate of some- 
thing less than 4. 

It is also interesting to note that some of the higher members of the 
paraffin series, such as n-heptane and in particular isooctane 

(CsHig), behave in solutions with perf luorocarbons as if their solubility 

parameters are significantly greater than those calculated from presumably 

V 2 

reliable AH density data at solution temperature. 

The effect is greatest with i-Cg H 18> which has a labile tertiary hydrogen 
atom that may behave somewhat as a Lewis acid. However, the effect apparently 
disappears with paraffins below 7 carbon atoms whose hydrogen atoms are less 
labile (at least toward thermal decomposition). Thus, such an effect would 
not be expected with methane. 

Negative deviations from ideal solution behavior would be expected pri- 
marily with strong interactions (e.g., dipole, acid-base, or hydrogen bonding) 
between solute and solvent, which appear to be unlikely to occur with methane, 
at least with common solvents. 

On the other hand, physical adsorption of methane on surfaces such as 
activated carbon or silica gel do in fact represent strong (van der Waals 
forces) interactions between methane and a chemical substrate (solid, in this 
case) which are comparable to negative deviations from "ideal interaction.” 
Thus, negative interactions with methane are possible and might be simulated 
in the liquid phase by surface forces either in mixed interacting solvent 
combinations or by adding solid surfaces to form a colloidal suspension. 

Data Collection 

The solubility data found in this study, covering a broad spectrum of 
solubility parameters and molal volumes as well as polarities, are shown in 
Table 2, and in Figure 1 plotted as a function of solubility parameter. In 
general, the data correlates reasonably well with the solubility parameter 
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Figure E-l. METHANE SOLUBILITY - SOLUBILITY PARAMETER CORRELATION 
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from 5.5 to 10 with solubility decreasing with increasing solubility param- 
eter® 


Most of the solubility data available in the literature were measured in 
terms of mole fraction of methane in equilibrium solution at the standard con- 
ditions of 25°C and 1 atmosphere total pressure. For our purposes, these data 
have been converted to molal solubilities at more practicable conditions of 34 
atmospheres (500 psia) and 25°C (77°F) assuming that Henry f s Law is obeyed, 
which is reasonably accurate for methane up to at least 34 atm. Solubilities 
in the lower paraffinic hydrocarbons (C 2 -C 4 ), however, were derived from 
vapor-liquid equilibrium data for those systems in which both components have 
appreciable, or high concentrations in both phases in two ways: 1 ) equilibrium 
concentrations in liquid at a total pressure of 34 atm as measured, and 2) by 
correcting the measured mole fractions of methane in the liquid phase at 25°C 
to a methane partial pressure of 34 atm via Henry’s Law® 

Solubility @ 34 atm CH^, 25°G = 


Mole Fraction in Liquid Phase @ 25 Q C x 34 (7) 

Mole Fraction CH^ in Vapor Phase x Total Press X 

Such values are fictitious in the sense that they can be achieved only at pres- 
sures much greater than 34 atm or not at all in the case of ethane. However, 
they do represent the approximate solubilities which would be achieved if the 
solvent had the same solution properties as C 2 -C 4 hydrocarbons but were non- 
volatile. 

It will be noted from Figure 1 that the solubilities of methane in a wide 

variety of solvents other than alcohols correlate more or less linearly with 

2 

the solubility parameter of the solvent on a log molal solubility versus S 
basis, semiquant itatively in agreement with Equations 3 or 4. However, Figure 
1 is intended primarily to be Illustrative and not as a definitive 
correlation. Thus, the scatter in the data may result from failure to account 
properly for differences in molal volume of solute and solvent as well as 
inaccuracies in: 

• ' Solubility measurements 

• Estimation of solubility parameters 
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m Deviations from ideal gas behavior. 

There are several features of the correlation that need comment. 

1. Applicability of Solubility Parameter Values 

Although the correlation conforms, in general, to Equation 3 which has a 
theoretical basis, it should not be considered rigorous and quantitative. In 
particular, the assignment of quantitative values to the solubility parameter, 
S, is somewhat in question since values of heats of vaporization and molar 
volumes at solution temperature (25°) are sometimes difficult to find. Fur- 
thermore, as Hildebrand has pointed out, in some cases paraffin hydrocarbons 
in perf luorocarbons tend to behave as if the solubility parameter of one (or 
bath) of the components is shifted. Thus, there is a tendency in the 
literature to use alternate methods of estimating S, in particular, by de- 
fining the solubility parameter using empirical solubility determinations with 
solvents or solutes whose solubility parameters are known. 

The solubility parameter of methane at 23°C cannot be estimated by the 
method used in this study (based on AHy + V^), since 25°C is considerably 
above its critical temperatures (Tq = -82°C). Thus, the heat of vaporization 
at 25 °C is presumably zero and the net energy of vaporization (AH V - RT) may 
actually be negative (i.e., (0 - RT) ■ -592 cal/mole). However, in our case 
we have skirted the problem by correlating in terms of rather than (S^ - 
S 2 ) as suggested by Equations 3 or 4. In any case, the solubility parameter 
of methane corresponding to the values used for various hydrocarbon solvents 
is presumed to be less than those estimated for propane (S = 5.85) from actual 
AH V + V M data, and for ethane (S =4.05) estimated by extrapolation of A Hy 
from lower temperatures and by extrapolation of S values at lower temperatures 
from measured AHy + data. Thus, the solubility parameter of methane is 
presumably less than 4.0, comparable only to the perf luorocarbons , silicones, 
or lower hydrocarbons which have been found to be the best solvents. 

2. Effect of Molecular Weight 

It will be noted that with the paraffin hydrocarbon series, the correla- 
tion appears to hold up to about Cy above which the mole fraction solubility 
increases with increasing molecular weight. These deviations are qualita- 
tively if not quantitatively accountable in terms of the differences in molal 
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volumes which affect the third term in Equation 4. The somewhat high solubil- 
ities in perf luoroheptane and octaraethyl cyclotetrasiloxane can also possibly 
be explained in these terms* However, the physical explanation of this phe- 
nomenon may be that in the larger molecular weight solvents, the dissolved 
methane is interacting not with the molecule as a whole but with smaller seg- 
ments of the molecule. This physical picture is consistent with the fact that 
(n-C^Fg)^N which might be expected to behave as a spherical entity appears to 
conform to the solubility parameters correlation very well, whereas n- C 7 F 16 
and octamethyl cyclotetrasiloxane do not. This concept could be important in 
designing higher molecular weight solvents having good capacity on a gravimet- 
ric as well as a mole fraction basis. It is unfortunate that with the 
n-paraffin series of solvents, the increase in solubility on a mole fraction 
basis with increasing molecular wieght above n-heptane is not sufficient to 
overcome the loss in gravimetric capacity due to the molecular weight increase 
(compare solubilities in Table 2 on a mole fraction and lb/lb basis)* 

On the other hand, if high solvency of propane or particularly methane, 
ethane, and propane could be simulated in terms of a multiplicity of methyl-, 
ethyl-, or propyl-group branches on a surface (e.g., silica gel) or on an 
oligomer, the increase in molecular weight might be overcome by the increase 
in solvency in the side chains so that the solvent maintains the methane 
capacity on a gravimetric basis close to that of methane, ethane, or propane. 
Ethane at 3.6 MPa vapor pressure, for example, would dissolve methane at a 
ratio of 0.35 g CH^/g ethane. If this solvency can be maintained with the 
ethyl group branches on a low polymer, a gravimetric capacity exceeding that 
of adsorption could be achieved. 

3. Aliphatic Hydroxyl-Group Effect 

The solubilities of methane in aliphatic alcohols, although correlatable 
with solubility parameters appear to be greater at the same solubility param- 
eter value than found for the other solvents having no hydroxyl groups. This 
suggests that hydrogen bonding, if strong enough, does contribute to solvent 
power for methane independent of the solubility parameter. However, the 
contribution of the hydroxyl groups to solubility tends to be overwhelmed by 
their effect on solubility parameter. Nonetheless, the very existence of the 
effect suggests that it should be possible to utilize it to design systems 
with enhanced solubility for methane as discussed below. 
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a. Chemical Modification of Adsorbents 


One approach would be to enhance the adsorption capacity of solid adsor- 
bents by providing a multiplicity of hydroxyl (or other hydrogen bonding 
group, including amides) on the surface. Such adsorbent surfaces might be 
expected to exhibit a higher heat of adsorption, so that greater surface 
coverage would be achieved at lower pressures than for adsorbents with non- 
polar surfaces. This would suggest the coverage of high surface area 
adsorbents with the specific groups (organic or inorganic) having the highest 
affinity for methane presumably by hydrogen bonding. 

This concept is consistent with the results of Chuik, et al. ,^ which 
showed that reaction of the silanol groups of a silica gel with an amino 
organo silane compound reduced both the surface concentration of silanol 
groups and the adsorption capacity for methane and other hydrocarbons. How- 
ever, these results are not straightforward, since the treatment replaced the 
silanol groups with rather bulky 


OCH, 

I I 3 

-Si-0-Si-(CH_ )«NH 0 

I I 232 

och 3 

groups, which could' have prevented adsorption by bulk alone, 
b. Multicomponent Liquid Solvent System 

Design multicomponent liquid solvent systems in which the hydroxyl or 
other hydrogen bonding function is built-in in such a way as to enhance sol- 
vency. One approach might be a dual colloidal system in which a discrete 
second liquid or solid phase material provides a multiplicity of hydroxyl 
groups, interfaced with and dispersed in a continuous phase having a low solu- 
bility parameter. The easiest example would involve a dispersion of colloidal 
silica (millimicron diameter) in a low solubilty parameter solvent such as a 
hydrocarbon, perf luorocarbon or silicone. Such a material could, of course, 
be made either solid or liquid. In this case, the dispersing function would 
be provided by a partial coverage (by reaction) of the surface with low solu- 
bility parameter tails compatible with the matrix phase (silicones). 
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c. Dual Function Solvents 

Design single, dual function molecules or dual functional mixtures of 
solvent molecules containing both a low solubility parameter function and a 
hydrogen bonding or other affinity function in the form of a hydroxyl, amido 
or other groups. In particular, if the hydroxyl group effect on the solubil- 
ity parameter of the solution could be minimized by either: 

• Hindrance of the hydroxyl effect on S by surrounding it with low solu- 
bility tails in a single solvent molecule. 

• Neutralization of the effect of the hydroxyl group on solubility param- 
eter of the solution by hydrogen bond formation in a manner similar to 
the effect of ether formation on solubility parameter and solvent power 
of alcohols (e.g., S for diethyl ether = 7.5 compared to 12.7 for ethyl 
alcohol from which it is made). Similarly perhaps, the solubility param- 
eter of R-0***H-CH3 may be considerably less than ROH itself. If so, 
solvent raiStures containing a low solubility parameter solvent molecule 
as a matrix plus a strongly hydrogen bonding species (such as alcohol) in 
concentrations comparable to that of the expected solute concentrations 
could exhibit the advantages of both hydrogen bonding and low solubility 
parameter. 

d. .Hydrogen-Bonding Groups Other than Hydroxyl 

Find other organic or inorganic groups having even greater hydrogen bonding 

power than aliphatic hydroxyl groups. In particular, Wolfenden has reviewed 

and analyzed the "hydrophilicity" of various organic groups in terms of the 

water-to-vapor distribution coefficients for various aqueous solutions of 

12 

uncharged organic compounds containing these groups. The results suggest 
that the hydrogen bonding propensity of various groups increases in the fol- 
lowing order: Alkanes < - SH < - Cl < ethers and thioethers < esters < alde- 

hydes and ketones < nitriles < amines < alcohols < ^0 < acids, phosphotri- 
esters < Amides and diols < q-substituted amides < peptides < guanidines. 

This suggests that q-substituted amides and guanidines in which the 
q-substituent is a long chain hydrocarbon or perf luorocarbon dissolved in an 
appropriate solvent of low solubility parameter should be worthy of evalu- 
ation. 

14 

Solubility of Methane in Polymers 

Considerable literature exists on the solubility of various gases including 
methane in a number of polymers. Most of this work was aimed at determining 
the diffusion characteristics of polymeric films, one of the basic parameters 
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of which is solubility* The best polymeric solvent found so far is silicone 
rubber, which exhibits a significantly greater capacity for methane (0,018 
g/g) than polyethylenes, polyisobutylene and other systems. 1 ^ Unfortunately, 
the gravimetric capacity of even silicone rubber is very low compared to the 
low molecular weight solvents. 

The factors controlling solubility of fixed gases in these systems was 
not completely delineated by the literature review. For example, the 
solubility of methane in polyethylene in some studies appears to be limited to 
(or correlatable with) the amorphous phase . Yet other studies have found 
relatively little difference between the highly crystallizable linear poly- 
ethylene and amorphous systems such as polyisobutylene. ^ The theory seems to 
be trending toward a "hole" theory comparable to adsorption in such common 
adsorbents as activated carbon in that the solubility is a function of the 
"hole" volume. ^ 

Nonetheless, solubility would be expected to be greatest in relatively 
low solubility parameter polymers such as polyisobutylene and the silicones. 
Furthermore, amorphous polymers would be expected to exhibit greater "hole 
volume" than crystalline polymers. In any case, one wonders if the free 
"hole" volume in polyisobutylene or silicone polymers could be controlled by 
controlled cross-linking while in a swollen state with an easily volatilized 
solvent. Of particular interest would be silicon-containing polymers in- 
cluding: 

• Polyvinyl trimethyl silane which has been shown to dissolve methane to 

the extent of four times that of polyethylene terephthalate. 1 ^ 

• Poly l-(Trimethyl Silyl)-l-propyne which has recently been shown to 

dissolve oxygen to an extent ten times that of other known polymers. 21 

It is interesting to speculate as well that specialized adsorbents, as 
opposed to absorbents, could be developed from polymeric systems cross-linked 
in the swollen state, containing hydrogen bonding groups such as hydroxyl or 
amide groups, e.g., polyvinyl alcohol or polyethylene terephthalate. Such 
polymeric systems are sometimes called reticulated polymers and offer the pos- 
sibility of adsorbents with surface chemistries tailored for the particular 
adsorbate desired. We do not have surface area data on such systems but dry 
reticulated ion exchange resins have been evaluated as catalysts. 
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Such systems could be of considerable interest, at least conceptually, 
since they may lead to negative deviations from ideal interaction at the solid 
surfaces. However, strong adsorption would be expected to occur only for 
raonomolecular surface coverage. Thus, surface area would be important. 
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SUMMARY AND CONCLUSIONS 


The available literature data on the solubility of methane in a wide var- 
iety of solvents has been collected and analyzed and the general conclusions 
are as follows: 

• In general, the data pretty well conforms to the classical views of 
Hildebrand on regular solutions, i.e., solutions "in which orienting and 
chemical effects are absent and in which the distributions and orienta- 
tions are random..." In other words, methane behaves in most solutions 
as non-polar molecules (including the rare gases) subject primarily to 
Van der Waals dispersion forces. 

• The primary parameters controlling solubility in regular solution theory 
are the molal volume and the Hildebrand solubility parameter relative to 
that of methane, the solubility parameter (or cohesive energy density) 
being defined as: 

s = (— | ) 1/2 

V 1 


where: 

AHj = Heat of vaporization at solution temperature 
= Molal volume of Component 1. 

• Since the solubility parameter of methane is at the lower end of the 
scale (S = <4), the best solvents for methane are the perfluorocarbons, 
the silicones, and the lower aliphatic hydrocarbons such as propane. 

• The best solvent for methane so far identified is propane, in which 
methane is soluble to the extent of 0.15 mole fraction or 0.063 g CH^/g 
solvent. However, this is a lower capacity on a mass ratio basis than 
can be expected from adsorption on activated carbon. 

• On a mole fraction basis, the solubility of methane in octamethyl cyclo- 
tetrasiloxane (0.32) and in perfluoro n-heptane (0.28) is greater than in 
propane but on a gravimetric basis the solubilities are much lower (0.025 
and 0.016 respectively) due to their high molecular weight. Thus, methane 
solubilities in lower molecular weight perfluorocarbons and silicones 
should be determined. 

• In the aliphatic hydrocarbon series, the solubility of methane decreases 
with increasing molecular weight (and increasing S) up to about CL 
(n-hexane). Above (L, however, the solubility of methane on a mole frac- 
tion basis in this class of solvents increases rapidly with increasing 
molecular weight to 0.32 for C 32 (squalane). Unfortunately, the solvent 
weight increases more rapidly than solubility so that the solubility on a 
gravimetric basis continues to decrease with molecular weight. 
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9 The solubility of methane in isobutane (0.061 g/g) and neopentane (0.057 
g/g) suggests that the solubility on a gravimetric basis may be main- 
tained while molecular weight increases for highly branched paraffin 
hydrocarbon. Unfortunately, data on the solubility of methane in higher 
molecular weight highly branched paraffins is lacking. 

• In particular, preliminary calculations suggest that it may be possible 
to design higher molecular weight solvents having gravimetric capacities 
in the range of, or greater than, for adsorption. The solubility data 
presented above for the lower paraffinic hydrocarbons (e.g., propane) 
actually represent the concentration of methane in the liquid phase of a 
mixture in equilibrium with a vapor phase containing a considerable 
proportion of solvent vapor. For example, if the equilibrium 
concentration of methane in propane at 25°C and 3.5 MPa total pressure is 
corrected to 3.5 MPa partial pressure of methane, the solubility would 
increase from 0.063 g/g to 0.11 g/g propane. With ethane, the equilib- 
rium concentration of methane in the liquid phase at 25°C and 4.9 MPa is 
about 0.05. Thus, if corrected to 3.5 MPa partial pressure of methane, 
the solubility in methane becomes 0.35 g/g of ethane. Thus, if we could 
tie down the methyl, ethyl or propyl groups and still achieve their sol- 
vent power for methane, it may be possible to achieve a viable non- 
volatile solvent for our purposes. In particular, if we could simulate 
the high solvency of methane, ethane, or propane in terms of a multiplic- 
ity of methyl, ethyl, or propyl groups as branches on a surface (e.g., 
silica gel) or an oligomer, the gravimetric solution capacity for methane 
might be maintained at a reasonable level. For example, if the solvency 
of ethane can be maintained with ethyl groups on a low molecular weight 
polymer, a gravimetric capacity exceeding that for adsorption might be 
achieved: 
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9 The solubilities of methane in aliphatic alcohols, although correlatable 
with solubility parameters, appear to be greater at the same solubility 
parameter value than for other solvents having no hydroxyl groups. This 
suggests that hydrogen bonding power can contribute to solvent power for 
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methane independent of the solubility parameter. However, this contribu- 
tion to solvency by the hydroxyl group is overwhelmed by its effect on 
solubilty parameter. Nonetheless, the effect suggests that it may be 
possible to: 

1) Enhance the adsorptive capacity of solid adsorbents by providing a multi- 
plicity of hydroxyl groups on the surface. 

2) Design multicomponent liquid solvent systems in which the hydroxyl or 
other hydrogen bonding function is build into the system to enhance 
solvency. 

3) Design dual function systems, either single detergent-like molecules or a 
dispersion of hydroxyl group containing molecules in a matrix of low sol- 
ubilty parameter, in which the hydroxyl group effect on the solubility 
parameter of the solvent as a whole is suppressed or minimized by. 

4) Evaluate the effect of other groups having stronger hydrogen bonding 
power than aliphatic hydroxyl groups including amides, n-substituted 
amides and n-substituted guanidines. 

5) Of particular interest would be tailor-made polymers (or oligomers) of 
various silicon containing monomers and particularly the trialkyl silanes 
(as opposed to the dialkyl silicone type). 
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Table 2. EXPERIMENTAL SOLUBILITIES OF METHANE 
IN VARIOUS SOLVENTS 


Solvent 


Ethane 

Propane 

n-Butane 

i-Butane 

n-Pentane 

i-Pentane 

neo-Pentane 

n-Hexane 

z-Methylperitane 

n-Heptane 

n-Octane 

224-Trimethl 

Pentane 

n-Decane 

Polyethylene 


Cyclohexane 

CIS 1,2 Dimethyl 

Cyclohexane 

Trans-Cyclohexane 

1,3 Dimethyl 

Cyclohexame 

CIS 1,4 Cyclohexane 

Methyl Cyclohexane 

Bicyclohexyl 

Cyclooctane 


Benzene 

Toluene 

0- Xylene 
m-Xylene 
p-Xylene 
Diphynyl 

Methane 

1- Methyl 

Naphthalene 


Ref. 

S 


ALIPHATIC 

38 

4.05 

1,2,3 

5.85 

4,5 

6.6 

6 

6.2 

7 

7.0 


6.8 

14 

6.2 

8-11,13,35 

7.3 

12 

— 

35 

7.45 

17 

7.65 

15 

6.85 

17 

7.7 

41 

— 


ALICYCLIC 

9 

8.2 

24 


24 


24 


24 


25 

7.8 

16 

8.2 

28 

8.5 


AROMATIC 

9,23 

9.15 

19,25,35 

8.9 

26 

9.0 

26 

8.8 

26 

8.8 

16 

9.1 


31,32 


Methane Solubility 
at 25°C 

34 atm CH/, Pressure, 
@ Mol 

Faction, g CH 4 /g 
CH^ Solvent 

HYDROCARBONS 


0.40 

0.21 

0.235 

0.11 

0.192 

0.66 

0.205 

0.071 

0.175 

0.175 

0.177 

0 . 

0.218 

0.062 

0.17 

0.038 

0.17 

0.038 

0.16 

0.030 

0.17 

0.029 

0.184 

0.032 

0.20 

0.028 

— — 

0.005 

HYDROCARBONS 

0.11 

0.024 

0.135 


0.145 


0.145 


0.147 


0.135 


0.115 


0.094 


HYDROCARBONS 

0.071 

0.012 

0.076 

0.013 

0.085 

0.014 

0.092 


0.098 

0.016 

0.061 


0.053 



Equilibrium 


Liquid 

Composition 


at 25°C 

and 34 

atm Total 

Pressure 

Weight 

Percent 

g/g 

0.05 

0.027 

0.148 

0.063 

0.17 

0.054 

0.18 

0.061 

0.046 

0.046 

0.165 

— 

0.206 

0.057 
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Solvent 


(C 4 F 9 ) n 
C 7 F 16 

Hexafluoro 

Benzene 

Octamethylcyclo 
Tetrasilane 
Silicone Rubber 


CHoCl-CHoCl 

ccl 4 

Chlorobenzene 

CHCI 3 

cci 2 f-ccif 2 


h 2 o 

Methyl Alcohol 
Ethyl Alcohol 
Propyl Alcohol 
i-Propyl Alcohol 
n-Butyl Alcohol 
n-Amyl Alcohol 
n-Hexyl Alcohol 
n-Heptyl Alcohol 
n-Octyl Alcohol 
n-Decyl Alcohol 
Cyclo Hexyl Alcohol 


Table 2. EXPERIMENTAL SOLUBILITIES OF METHANE 
IN VARIOUS SOLVENTS, CON'T. 


Methane Solubility 
at 25°C 

34 atm CH^ Pressure, 
@ Mol 

Faction, g CH 4 /g 
Ref . S CH^ Solvent 

PERFLUOROCARBONS AND SILICONES 


21 

5.9 

0.23 


22 

5.7 

0.28 

0.016 

23 

8.25 

0.13 


31 

5.5 

0.32 

0.025 

42 



0.018 


CHLORINATED 

HYDROCARBONS 

35 

8.6 

0.029 

0.097 

0.011 

9,35 

9.5 

0.068 


35 

9.25 

0.058 

0.008 

15 


0.17 



Equilibrium 
Liquid 
Composition 
at 25 °C and 34 
atm Total Pressure 
Weight 

Percent g/g 


ALIPHATIC ALCOHOLS 


23.4 

0.0009 

14.5 

0.029 

12.7 

0.044 

11.9 

0.055 

11.5 

0.049 

11.0 

0.065 

10.3 

0.073 

10.0 

0.080 

9.5 

0.088 

9.3 

0.095 

0.108 

0.043 
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Table 2. EXPERIMENTAL SOLUBILITIES OF METHANE 
IN VARIOUS SOLVENTS, CON'T. 


Methane Solubility 
at 25°C 

34 atm CH^ Pressure, 
@ Mol 

Faction, g CH^/g 

Solvent Ref. S CH^ Solvent 



MISCELLANEOUS OXYGEN 

COMPOUNDS 

Diethyl Ether 

9 

7.45 

0.154 ( 

Dioxane 

9 

10.0 

0.044 

Methyl Acetate 

9 

9.5 

0.066 

Acetone 

9 

9.65 

0.063 

Ethylene Oxide 

27 


0.055 

Methyl Ethyl 
Ketone 


9.3 

0.06 

Propylene 

Carbonate 

39 

13.3 

0.023 

Glycerol tri- 
Acetate 

39 


0.065 

Sulfur Dioxide 
COS 

29 


0.035 

0.076 

cs 2 

21,22 

10.0 

0.044 

Dimethyl 

Sulfoxide 

30 


0.013 

Tributyl 

Phosphate 

39 


0.155 


MISCELLANEOUS NITROGEN COMPOUNDS 

Ammonia 

33 


0.012 

Aniline 

35 

11.6 

0.023 

Cyclohexyl 

Amine 

29 

8.7 

0.065 

Methyl 

Pyrrolidane 

39 


0.033 

Nitrobenzene 

35 

10.0 

0.030 
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